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I  SUMMARY 


This  report  describes  the  work  done  on  the  development 
of  a  deep  hole  stress  measurement  device  under  Contract  No. 
R44620-7 2-C-0029 .  Our  intention  was  to  field  test  the  device, 
however,  we  did  not  get  that  far.  In  this  report  we  describe 
the  mechanical  design,  the  development  of  the  telemetrv  and 
control  systems,  and  the  considerable  effort  spent  on  studying 
the  problems  associated  with  residual  stresses  in  in  situ 
rocks  and  their  effects  on  stress  measurement  techniques. 

Data  reduction  from  the  device  has  also  received  a  gfeat 
deal  of  attention.  Three  different  approaches  were  tried  and 
these  are  also  described  in  the  report. 
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II .  Introduction 

The  purpose  of  this  investigation  was  to  develop,  build 
and  test  an  instrument  to  accurately  measure  primary  rock 
stresses  in  the  side  walls  of  long  (1000  meters  or  more)  6*5 
inch  diameter  drill  holes.  This  instrument  would  be  used  to 
measure  stresses  around  existing  uncased  oil  field  holes. 

In  particular  we  wanted  to  be  able  to  make  this  instrument 
available  for  projects  such  as  the  Rangely,  Colorado,  Oil 
Field  Earthquake  Investigation.  In  this  and  other  earth¬ 
quake  investigations  one  of  the  most  important  missing  bits 
of  information  required  for  complete  analysis  of  the  problem 
is  the  existing  primary  state  of  stress  in  the  ground.  Other 
potential  areas  of  application  for  our  instrument  include: 


site  investigations  for  nuclear  power  plants,  determining 
suitability  of  a  site  for  underground  nuclear  testing,  in¬ 
vestigating  an  underground  location  as  a  potential  storage 
cavern  for  L.  N.  G.  (liquefied  natural  gas),  site  investi¬ 
gations  for  major  tunnels  and  underground  power  houses,  and 
planning  studies  for  deep  underground  mines.  At  the  present 
we  know  of  no  existing  technique  or  apparatus  which  can 
accurately  determine  the  three  principal  stresses  that  are 
acting  and  their  orientations  in  long  drill  holes. 

State  of  the  Art 

Numerous  rock  stress  measurement  techniques  have  been 
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developed  based  on  the  principle  of  strain  relief.  The 
general  idea  is  that  the  rock  is  believed  to  be  under  some 
initial  stress.  Some  of  the  rock  is  removed  and  deformation 
of  the  remaining  rock  occurs.  The  deformation  is  measured 
and  recorded.  Usually  the  next  step  is  to  make  the  assump¬ 
tion  that  the  rock  is  isotropic  and  linearly  elastic.  If 
the  geometry  of  the  strain  relief  system  is  simple  enough, 
an  analytical  solution  based  on  the  assumptions  of  linear 
isotropic  elasticity  is  found.  This  sol  ition  relates  primary 
or  field  stresses  to  the  strains  in  the  rock  where  the 
deformation  was  measured  through  the  elastic  moduli  of  the 
rock.  The  moduli  arc  determined  by  laboratory  tests  on 
samples  of  the  rock  and  these  moduli  along  with  the  measured 
strains  and/or  deformations  arc  put  into  the  equations  of  the 
analytical  solution  and  the  equations  are  solved  for  the 
primary  stresses. 

If  the  geometry  of  the  system  is  too  complicated  for  an 
analytical  solution  to  be  found  there  arc  still  three  methods 
of  getting  the  primary  stress-strain  at  the  measuring 
position  relationship:  (1)  a  numerical  approximation  can 
be  made  using  finite  difference  of  finite  element  techniques. 
This  procedure  is  usually  limited  in  practice  to  two- 
dimensional  analyses  because  of  limitations  on  the  size  of 


the  available  computers;  (2)  a  three-dimensional  frozen 
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stress  photoelasticity  model  is  made  and  analyzed.  The 
major  criticism  of  this  method  is  that  Poisson's  ratio  in 
the  model  material  is  0.50  whilst  in  rock  it  is  commonly 
around  0.20;  and  (3)  full  scale  laboratory  calibration  tests 
can  be  made  in  large  blocks  of  rock  with  known  stresses 
applied  by  a  loading  frame. 

Several  different  techniques  have  evolved  in  the  past 
20  years  for  making  strain  relief  stress  measurements.  These 
techniques  differ  from  each  other  only  in  the  details  of 

(1)  how  they  make  the  deformation  or  strain  measurement, 

(2)  how  the  measuring  components  are  arranged  relative  to 
the  strain  relieving  activity,  and  (3)  just  how  the  strain 
relief  is  accomplished.  Two  of  the  most  widely  used  tech¬ 
niques  are  the  borehole  deformation  gage  and  the  borehole- 
end  strain  relief  methods  or  "doorstoppers." 

Both  borehole  deformation  gages  and  "doorstoppers"  have 
been  extensively  tested  in  the  laboratory  and  successfully 
used  in  the  field.  There  are  many  advantages  and  disadvan¬ 
tages  claimed  for  each  technique  but  the  main  problem  with 
either  of  them  is  that  the  maximum  depth  of  hole  that  they 
can  practically  be  used  in  is  of  the  order  of  50  to  100  feet. 
Another  disadvantage  common  to  both  techniques  is  that  a 
drilling  rig  must  be  set  up  on  the  site  and  special  holes 
drilled  for  the  in  situ  stress  determinations. 
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Another  candidate  technique  for  deep  hole  in  situ 
stress  measurements  is  hydraulic  fracturing.  Hydraulic 
fracturing  is  a  technique  originally  developed  to  stimulate 
production  from  oil  wells.  A  section  of  a  borehole  is  sealed 
with  straddle  packers  and  fluid  pressure  is  then  applied  to 
the  bare  walls  of  the  hole  between  the  packers.  The  pressure 
is  increased  until  a  fracture  develops  in  the  wa] 1  and  the 
fracture  propagates  outwards  as  the  fluid  flows  into  it. 
Propagation  is  perpendicular  to  the  direction  of  the  minimum 
principal  stress,  and  the  fluid  pressure  required  to  extend 
the  fracture  is  approximately  equal  to  the  minimum  principal 
effective  stress  in  the  rock.  This  direct  measure  of  the 
minimum  principal  stress  is  independent  of  rock  properties. 

Kehlo  (1964)  has  given  a  more  refined  analysis  to 
determine  the  fluid  pressure  required  to  initiate  a  fracture 
in  rock  around  a  circular  borehole.  This  analysis  is  based 
on  the  usual  assumptions  that  the  rock  is  elastic,  isotropic, 
and  homogeneous.  In  addition  the  rock  must  not  be  jointed 
or  faulted.  Fairhurst  (1968)  has  used  this  analysis  to 
attempt  to  determine  all  three  principal  stresses  and  their 
directions  from  carefully  conducted  hydraulic  fracture  tests 
and  subsequent  observations  of  the  orientations  of  the 
fractures  that  wore  formed.  Hydraulic  fracturing  can  be 
done  in  holes  of  any  depth,  and  it  can  be  done  in  any 


existing  open  (uncased)  hole.  The  principal  limitation  of 
the  technique  is  that  it  is  difficult  in  practice  to  get  much 
more  than  a  single  (the  minimum)  principal  stress  magnitude 
from  the  field  measurements.  Haimson  (1973)  and  (1974) 
and  Healy,  et.  ai  .  (1974)  have  given  recent  examples  of 
stress  determination  by  hydraulic  fracturing  techniques. 

The  technique  which  we  have  attempted  to  develop  in 
this  project  is  based  on  the  strain  relief  principal  but  is 
one  which  should  be  able  to  be  used  in  existing  uncased 
holes  of  any  depth.  We  have  run  into  several  unexpected 
problems  in  the  course  of  this  research  and  are  not  nearly 
as  far  along  as  we  had  proposed  or  intended  to  be  at  this 
time.  We  have  not  yet  run  into  insurmountable  problems, 
however,  and  we  do  intend  to  continue  working  on  the  tech¬ 
nique  until  it  is  a  reliable  field  stress  logging  tool. 

General  Description  of  the  Apparatus 

We  have  designed,  built  and  laboratory  tested  a  device 
capable  of  determining  the  complete  state  of  stress  in  the 
ground  from  measurements  made  in  the  side  walls  of  6  to  7 
inch  diameter,  uncased  oil  field  holes.  The  principle  we 
use  is  strain  relief  undercoring  using  friction  bonded  strain 
gages.  This  technique  was  first  proposed  by  Hoskins  (1968) 
for  use  at  the  end  of  long  but  smaller  ( 2  to  3  inch)  diameter 


boreholes.  In  a  6  to  7  inch  hole  there  is  enough  room  to  make 


the  measurements  in  the  borehole  walls  and  this  is  a  con¬ 
siderable  advantage.  No  borehole  end  preparation  is  required 
and  stresses  can  be  determined  at  various  depths  in  the 
hole  by  repeated  tests. 

The  principle  of  strain  relief  stress  measurements  by 
drilling  a  small  hole  in  the  center  of  a  strain  gage 
rosette,  or  trepanning,  is  well  established  and  has  been 
widely  used  in  experimental  stress  analysis  (c.f.  Hetenyi, 
Handbook  of  Experimental  Stress  Analysis  1950) .  Hoskins 
(1968)  has  previously  experimentally  determined  stress 
concentration  factors  for  interpreting  results  of  this  type 
of  test  at  the  flattened  end  of  a  borehole.  The  same  sort 
of  experiments  were  done  as  a  part  of  this  project  for  tests 
performed  in  the  walls  of  a  borehole.  Since  an  elasticity 
solution  exists  for  the  stresses  in  the  walls  of  a  cylin¬ 
drical  hole  in  a  general  three-dimensional  stress  field, 
only  the  effects  of  the  strain  gage  rosette  and  trepanning 
drill  geometry  have  had  to  be  determined.  Careful  but 
straight-forward  laboratory  experiments  were  required. 

Based  on  the  requirements  outlined  in  the  introduction 
and  the  principal  investigator's  previous  experience  with 
other  in  situ  stress  and  strain  measuring  devices  we  adopted 
the  following  set  of  specifications. 
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Resolution:  10  ppm  strain  -  equivalent  to  approxi¬ 

mately  3  bars  in  rock  with  a  modulus  of 
5 

3  x  10  bars 

Operating  Depth:  2500  meters 
Hole  Diameter:  inches  +  h  inch 

Maximum  Temperature:  100°C 
Case  Material:  Stainless  Steel  (416) 

Orientation:  Sperry  Sun  Magnetic  Device 

Case  Size:  5k  inches  diameter 

Hydraulic  Power:  Bled  from  self-contained  intensifier 
Electrical  Power:  From  self-contained  battery  pack 
A  general  schematic  of  the  apparatus  is  shown  in 
Figure  1. 

Ill  .  Mechanical  Design 

A  brief  description  of  the  apparatus  follows.  There 
are  three  strain  rosette  trepanning  units.  They  are  mounted 
on  hydraulic  cylinders  120  apart  about  the  longitudinal  axis 
of  the  device.  There  are  si  .  other  hydraulic  shoes,  three 
above  and  three  below  the  measurement  units  to  firmly  fix 
the  entire  device  in  the  hole.  The  hydraulic  pressure  is 
generated  in  the  downhole  apparatus  and  controlled  from  the 
surface.  Rotation  and  advance  of  the  trepanning  drills  is 
done  by  small  electric  motors  also  contained  in  the  device 
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but  controlled  from  the  surface.  The  entire  apparatus  is 

constructed  as  a  waterproof  bomb  so  that  the  electric 

motors  and  the  data  transmission  package  are  operating  in 

a  dry  environment.  The  strain  gage  signals  will  be  fed  up 
a  cable  to  the  surface  and  recorded  on  standard  equipment. 
Power  to  the  drilling  motors  and  hydraulic  pumps  comes  from 
a  self-contained  battery  pack.  A  standard  borehole  survey¬ 
ing  instrument  will  be  fixed  to  the  stress  measurement  device 
to  indicate  the  orientation  of  the  strain  gages  relative  to 
geographical  coordinates.  A  major  feature  of  this  system 
is  that  no  drilling  rig  or  derrick  is  required  over  the 
hole,  only  a  wire  line  truck  carrying  the  cable  reel,  power 
supplies,  and  read  out  instrumentation.  We  will  not  have  to 
try  to  transfer  signals  up  or  down  a  rotating  drill  string. 

The  main  portion  of  the  apparatus  was  machined  from  416 
stainless  steel  bar  stock.  Design  emphasis  was  placed  on 
simplicity  for  the  machine  shop  and  extreme  overall  rugged¬ 
ness  of  the  final  package. 

Friction  Bonded  Strain  Gages 

The  friction  bonded  strain  gages  are  the  most  crucial 
element  in  the  entire  device.  It  was  necessary  for  us  to 
design  and  fabricate  our  own  gages  because  the  commercially 
available  gages  (Model  CBF-6  manufactured  by  Tokyo  Sokki 
Kenkyujo  Col,  Ltd.)  could  not  be  bent  to  conform  to  the 
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curvature  of  the  sides  of  the  borehole  or  adequat  ely  water¬ 
proofed.  Several  different  versions  of  the  friction  bonded 
gages  were  constructed  and  tested.  The  best  so  far  in  terms 
of  both  performance  and  ease  of  construction  is  a  composite 
consisting  of  3  nylon  plugs  3/4  of  an  inch  diameter  set  into 
a  steel  holder  with  a  backing  plate.  The  strain  gages  are 
bonded  to  the  nylon  plugs  with  their  leads  fed  through  the 
backing  plate  into  the  main  body  of  the  apparatus.  The 
outer  surface  of  the  gages  is  coated  with  a  thin  layer  of 
60  grit  carborundum  powder  in  an  epoxy  matrix.  This  outer 
surface  will  be  in  contact  with  the  side  walls  of  the  bore¬ 
hole  . 
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Telemetry  and  Control  System 

Telemetry  System 

The  telemetry  package,  except  for  strain  gauges  and 
instrumentation  amplifiers  was  supplied  by  IED,  a  division 
of  Conic  Corporation,  San  Diego,  California.  The  system 
is  capable  of  accepting  16  channels  of  analog  data  and 
converting  them  to  digital  form  for  transmissions  to  the 
surface  on  two  wires  of  the  seven  conductor  cable.  At 
the  surface  the  decoder  drives  a  digital  printer. 

Control  System 

The  control  system  is  a  digital  system  actuated  by 
a  5-position  Command  Selector  Switch.  The  commands  are: 
Telemetry  power,  Jacks  Out/In,  Strain  gauge  Jacks  Out/ 

In,  Drill  power,  and  Drill  Out/in.  The  command  signals 
are  encoded  in  binary  form  and  sent  over  4  wires  of  the 
•cable  to  the  downhole  decoder  and  memory  unit.  This  unit 
contains  a  decoder  and  various  relays  which  supply  power 
to  the  appropriate  control  valves  and  to  the  electric 
drill  and  telemetry  system.  These  power  signals  are  also 
sent  to  the  Status  Signal  Generator  which  sends  a  suit¬ 
able  status  signal  to  a  status  light  on  the  surface. 

A  pressure  transducer  is  used  to  sense  pressure  variations 
resulting  from  the  operation  of  control  valves  for  Jacks 
Out,  Strain  Guage  Jacks  Out,  and  Drill  Out. 


-----  ■  '  •  - 
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IV  Residual  Stress  Problems 

Any  in  situ  stress  measurement  technique  actually 
measures  the  sums  of  several  independent  stress  fields 
that  are  present .  These  independent  stress  components 
include:  (1)  the  stress  component  due  to  gravity,  (2) 
stresses  due  to  thermal  gradients,  (3)  the  stress 
component  due  to  the  penetration  of  the  rest  of  the  stress 
field  by  the  emplacement  hole  and  measurement  device,  (4) 
the  stress  component  due  to  currently  active  tectonic 
processes,  and  (5)  the  stress  component  due  to  locked 
in  or  residual  stresses.  In  order  to  make  intelligent 
use  of  the  in  situ  measurements  we  have  to  be  able  to 
separate  these  various  stress  fields  and  report  them 
independently.  The  vertical  component  of  the  gravity  stress 
is  known  as  well  as  the  density  of  the  rock,  the  vertical 
depth  below  a  free  surface,  and  the  gravitational  constant 
are  known.  It  can  be  calculated  from  equation  (1): 

(i) 

Where  <JZ  is  the  vertical  stress  due  to  gravity,  /> 

is  rock  density,  g  is  the  gravitational  constant  and  h  is 
the  depth  of  the  measurement  below  the  free  surface.  The 
horizontal  component  of  the  gravity  stress  field  depends 
upon  material  properties  and  material  behavior  as  well. 

If  the  rock  is  considered  to  be  homogeneous,  isotropic 
and  linearly  elastic,  the  horizontal  components  of  the 
gravity  stress  field  are  equal  and  can  be  calculated  by 
equation  (2)  : 

(2) 

Where  cr"H  is  the  horizontal  stress  due  to  gravity,  or  is 


Poisson's  ratio  and  C Tv  is  the  vertical  stress  due  to 
gravity.  Since  rocks  commonly  have  Poisson's  ratios  of 
0.20  to  0.25  this  equation  leads  to  the  usual  estimate  of 
the  horizontal  stresses  being  equal  to  1/4  to  1/3  of  the 
vertical  stress. 

Stresses  due  to  thermal  gradients  can  likewise  be  cal¬ 
culated  from  conventional  thermo-elastic  relationships  if 
the  magnitude  of  the  gradients  and  the  appropriate  physical 
constants  are  known. 

The  laboratory  and  theoretical  studies  normally  per¬ 
formed  to  prove  any  stress  measurement  device  before  it  is 
taken  into  the  field  include  determination  of  the  secondary 
stress  field  due  to  the  measurement  technique  itself. 

We  are  left  then  with  tectonic  and  residual  stresses 
to  determine  separately.  So  far  as  we  know  there  are  as 
yet  no  published  results  of  in  situ  measurements  which 
accomplish  this  separation. 

The  tectonic  stress  field  cannot  be  accurately  cal¬ 
culated  at  our  present  stage  of  urtferstanding  of  tectonic 
processes.  In  fact  it  is  the  tectonic  component  that  we 
are  usually  attempting  to  determine  when  we  make  in  situ 
measurements . 

The  residual  stress  field  is  here  defined  as  the  stress 
field  remaining  in  a  specimen  in  the  absence  ox  thermal 
gradients  in  the  specimen  or  external  loads  applied  to  its 
boundaries.  This  is  a  standard  definition  of  residual  stres 
as  used  in  metallurgical  and  ceramic  engineering  practice. 
Obviously  the  body  must  satisfy  internal  and  external  equi¬ 
librium  and  the  integral  of  stresses  taken  over  the  volume 
of  the  body  must  equal  zero. 
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While  residual  stresses  in  metals  have  lorn  een  known, 
routinely  measured,  and  manipulated  it  is  only  in  the  past 
few  years  that  they  have  received  any  detailed  attention  by 
workers  in  the  field  of  rock  mechanics.  Since  their  magni¬ 
tudes  can  approach  the  yield  strength  of  the  material  they 
cannot  be  dismissed  as  of  trivial  importance.  In  fact,  some 
of  the  high  in  situ  stress  values  reported  from  surface 
measurements  in  supposedly  stable  or  shield  areas  may  be 
mainly  residual  with  little  or  no  tectonic  component. 

Further,  since  the  blocks  of  rock  which  we  are  using  to 
laboratory  calibrate  our  deep  hole  device  do  in  fact  contain 
significant  residual  stresses  as  does  the  Weber  sandstone 
which  is  the  reservoir  rock  at  Rangely,  Colorado,  it  is  quite 
important  that  we  become  able  to  distinguish  residual  from 
tectonic  stresses.  Accordingly  we  have  devoted  considerable 
time  and  effort  to  the  study  and  measurement  of  residual 
stresses  in  various  rocks  in  the  laboratory  as  a  part  of  this 
Droj  ect . 

The  initial  experimental  work  was  reported  in  the  Final 
Report  on  Contract  No.  F44620-70-C-0073  entitled  "Development 
of  a  Deep  Hole  Stress  Measurement  Device".  Some  additional 
experimental  work,  a  discussion  of  the  mechanism  of  formation 
of  residual  stresses,  and  a  preliminary  attempt  at  numerically 
modeling  residual  stresses  have  been  performed  during  the 
present  contract.  The  work  was  reported  in  paper  by  J.E. 
Russell  and  E.R.  Hoskins  at  the  14th  Symposium  on  Rock 
Mechanics  and  the  following  sections  are  largely  taken  from 
that  paper. 

Fundamental  Concepts 

Perhaps  the  most  fundamental  concept  necessary  in  the 
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understanding  of  residual  stress  phenomena  is  that  of  the 
self-equilibrating  unit.  Since  residual  stresses  exist 
in  the  absence  of  loads  and  of  thermal  gradients,  they 
must  somehow  form  self-equilibrating  units  in  order  to  satisfy 
Newton's  second  law.  Self-equilibrating  units  that  are  either 
macro  or  micro  in  nature  are  possible  and  residual  stresses  are 
sometimes  classified  as  being  either  macro  or  microresidual 
stresses.  Examples  of  both  types  of  residual  stresses  are 
discussed  in  a  later  section. 

A  well  known  result  from  thermoelasticity  is  that  an 
increase  in  temperature  results  in  an  increase  n  volume. 

If  the  material  is  isotropic  and  homogeneous,  a  uniform 
change  in  temperature  will  cause  a  volume  change  with  no 
accompanying  change  in  the  stress  field  or  the  shape  of 
the  body.  On  the  other  hand,  temperature  gradients  induce 
both  changes  in  volume  and  in  shape.  If  the  change  in  volume 
‘or  shape  is  constrained  or  would  induce  a  displacement  field 
that  is  geometrically  incompatible,  a  thermally  induced  stress 
field  results  in  order  to  maintain  compatible  deformation.  The 
same  situation  exists  regardless  of  the  mechanism  that  induces 
the  inhomogeneous  strain  field.  Thio  concept  is  used  to 
explain  some  of  the  proposed  mechanisms  for  residual  stresses 
in  rock. 

The  terms  residual  stress  and  residual  strain  are  often 
used  interchangably  although  by  definition  they  are  not  the 
same.  The  term  residual  st  -ain  used  in  this  context  refers 
to  the  strain  field  maintained  by  the  residual  stress  field 
and  should  not  be  confused  with  the  permanent  set  (or  strain) 
induced  by  loading  a  material  beyond  its  yield  point  and  then 
releasing  the  loud  to  zero.  This  permanent  set  can  exist  (at 
least  ideally)  with  no  stress  while  residual  strain  cannot. 
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This  situation  is  shown  on  Fig.  2. 

The  presence  of  both  residual  stress  and  strain  fields 
leads  to  the  existence  of  internal  residual  strain  energy. 
Friedman  reports  th  t  the  internal  residual  strain  energy  is 
commonly  found  to  be  of  the  order  of  104  ergs/cm3.  These 
energies  arc  computed  us.tr  residual  strains  measured  in 
X-ray  if  fraction  studies  ^ 1  by  strain  relief  methods  along 
with  nc  appropriate  modulus  of  elasticity.  As  will  be  shown 
later  there  is  preliminary  evidence  that  at  least  part  of 
the  intern  1  residual  strain  energy  in  rocks  is  temperature 
dependent . 


Existence  of  Residual  Stresses  and  Strains  in  Rock 
Evidence  of  the  existence  of  residual  stresses  and  strains 
in  rock  exists  both  in  the  field  and  from  laboratory 
measurements.  In  the  field,  evidence  such  as  the  exfoliating 
of  practically  unweathered  rock  which  is  isolated  from 
active  tectonic  movements  has  been  reported  by  Varnes . 

Further  field  evidence  is  supplied  by  the  distortion  of 
saw  cuts  and  drill  holes  beyond  what  would  normally  occur 
elastically  under  conditions  where  the  only  reasonable 
explanation  is  locked  in  stresses. 

Laboratory  evidence  for  the  existence  of  residual 
stresses  and  strains  has  been  found  using  both  X-ray 
methods  and  strain  relief  methods.  In  addition,  optical 
evidence  such  as  undulatory  extinction  in  quartz  points 
toward  the  existence  of  residual  deformation  in  crystals. 

Mechanisms _ Producing  Residua]  Stresses  in  Rock 

Several  mechanisms  producing  residual  stresses  in  rocks 
have  been  discussed  in  the  references  given  in  this  section. 
These  mechanisms  are  briefly  reviewed  here  for  completeness 
and  to  point  out  their  relationship  to  the  fundamental 


concepts  discussed  above.  Some  variations  of  previously 
suggested  mechanisms  are  considered.  Most  of  the  mechanisms 
discussed  are  micro  rather  than  macro  in  nature.  It  should  be 
noted  that  the  designation  micro  as  used  here  refers  to  the 
size  of  the  self-equilibrating  unit  which  may  be  of  the  order 
of  several  times  the  average  grain  size. 

The  classification  of  a  mechanism  as  micro  may  be 
further  subdivided  into  micro-elastic  or  micro-dissipative 
depending  on  whether  elastic  elements  alone  or  elastic 
elements  in  conjunction  with  friction  elements  or  viscous 
elements  are  used  in  the  model  of  the  mechanism. 

Micro-elastic  residual  stresses  develop  when  a  rock 
containing  minerals  having  different  and/or  anisotropic 
coefficients  of  thermal  expansion  undergoes  a  uniform 
change  in  temperature.  Here  the  source  of  lionuniforin  strain 
is  the  nonuniformity  of  the  expansion  characteristics  coupled 
with  a  uniform  change  in  temperature.  This  mechanism  might 
be  referred  to  as  micro-thermoelastic .  It  should  be  noted  that 
since  the  residual  stresses  and  strains  in  this  case  vary 
with  temperature  the  internal  residual  strain  energy  will  also 
be  temperature  dependent.  Further,  this  mechanism  suggests 
that  strain  relief  measurements  made  on  an  unconnected  piece 
of  rock  at  depth  where  the  temperature  is  higher  will  probably 
be  different  than  strain  relief  measurements  made  on  the  same 
piece  of  rock  in  the  laboratory.  Some  preliminary  evidence 
pointing  to  the  existence  of  this  mechanism  in  a  medium 
grained  granite  is  reported  in  a  later  section. 

Two  other  examples  of  micro-elastic  residual  stresses 
in  rocks  are  provided  by  the  crystallization  of  qranitr 
under  pressure  and  the  cementation  of  sand  grains  under  high 
hydrostatic  pressure.  In  both  cases,  the  mechanism  is 


essentially  the  same  although  the  scale  /nay  be  different. 
Consider  the  very  simple  qualitative  model  shown  on  Fig. 

3.  Two  springs  having  different  spring  constants  are 
compressed  by  the  same  amount  and  are  initially  independent 
of  each  other.  If  we  now  connect  both  springs  to  the  same 
cross  bar  and  require  that  the  bar  can  move  only  parallel 
to  itself,  we  see  that  the  situation  is  unchanged  as  long  as 
the  loading  condition  remains  static.  However,  if  the  load 
is  removed  both  elastic  springs  attempt  to  regain  their 
respective  initial  unstrained  positions.  This  is  now 
impossible  because  of  the  newly  created  bond  between  the 
springs  and  tension  will  develop  in  the  stiffer  element  while 
the  more  flexible  element  will  remain  in  compression. 

On  a  smaller  scale,  the  self-balanced  stress  fields 
around  dislocations  in  the  crystal  lattice  arc  examples 
of  elastic  residual  stresses.  Dislocations  have  been  the 
subject  of  a  considerable  amount  of  research  and  much 
information  is  available.  On  the  smallest  scale,  microresidual 
stresses  are  due  to  misfitting  solute  atoms  and  individual 
dislocations.  c.f.  McClintock  and  Argon  (1966) 

Chemically  induced  volume  changes  can  be  the  source 
of  nonhomogencous  strain  that  induces  residual  stresses. 

This  mechanism  is  sometimes  mentioned  but  to  the  authors' 
knowledge  has  not  been  discussed  in  any  detail. 

Microrcsidual  stresses  may  also  involve  dissipative 
elements  such  as  friction  blocks  and/or  dashpots.  As  shown 
in  Fig.  4,  a  self-equilibrating  element  may  consist  of  a 
spring  and  a  friction  block  in  a  series  connected  to  a 
spring  in  parallel.  If  this  unit  is  loaded  beyond  the 
force  necessary  to  slide  the  friction  element  and  then 
subsequently  unloaded,  the  spring  in  scries  with  the  friction 
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Fig.  3:  A  purely  clastic  residual  stress  mechanism 
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Fig.  4:  An  elastoplastic  residual  stress  mechanism 
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element  will  develop  tension  while  the  spring  on  the  other 
side  will  remain  in  compression.  In  this  case,  the  residual 
strain  field  would  not  be  time  dependent.  On  the  other  hand, 
if  a  viscous  element  has  been  used  rather  than  a  friction 
element,  the  resulting  strain  recovered  would  be  time  dependent. 
Some  evidence  of  time  dependent  strain  relief  has  been 
reported  by  Varnes  (1969)and  Emery  (1964). 

On  the  macro  scale,  which  is  taken  here  to  be  orders 
of  magnitude  greater  than  grain  size,  it  is  difficult  to 
imagine  that  many  of  the  above  discussed  mechanisms  are 
operative  because  of  the  fragmented  nature  of  much  in  situ 
rock.  Nevertheless,  a  macro-dissipative  residual  stress 
mechanism  usually  assumed  to  operate  in  ductile  metals  may 
be  useful  in  explaining  some  phenomena.  Consider  a  plate 
flexed  by  some  external  agent  to  the  point  that  yielding 
occurs  on  the  upper  and  lower  surfaces  to  a  depth  less  than 
half  the  thickness  of  the  plate.  Then  an  elastic  core  remains 
around  the  neutral  surface  of  the  plate.  Upon  unloading,  the 
elastic  core  attempts  to  straighten  out  and  regain  its  original 
unflexed  position  which  is  now  impossible  because  of  the 
yielding  in  upper  and  lower  regions  of  the  plate.  This 
situation  is  shown  on  Fig.  5. 


It  should  be  pointed  out  that  it  is  likely  that  more 
than  one  of  the  above  mechanisms  may  he  acting  simultaneously 
in  any  given  rock.  For  example,  the  thermoelastic  mechanism 
may  be  acting  concurrently  with  nearly  any  of  the  others,  or 
macro  induced  microresidual  stresses  are  apparently  possible. 
The  possible  coexistence  of  several  mechanisms  complicates  the 
job  of  attempting  to  isolate  their  effects. 

Methods  for  Measuring  Residual  Strains 

Two  basic  techniques  used  for  measuring  residual  strains 


strain  Stress 


'•  Partially  yielded  plate. 
Tension 


S^Compressjory 
Partially  yielded  plate  unloaded 
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are  X-ray  diffraction  and  strain  relief  methods.  Measurement 
of  residual  strain  by  X— ray  diffraction  is  discussed  by 
Friedman,  (1968)  and  will  not  be  considered  here.  The  strain 
relief  methods  consist  of  removing  some  of  the  rock  by  sawing, 
coring,  etching  or  any  other  means  and  recording  the  chanqe 
in  strain  in  the  remaining  material.  Common  techniques 
include  overcoring  and  undercoring  (trepanning)  to  relive  the 
residual  stress  field.  The  amount  of  strain  relieved  is 
usually  recorded  by  bonded  electrical  resistance  strain  gages 
although  photoelastic  coatings  have  also  been  used.  The 
strain  changes  recorded  in  the  overcoring  case  may  be  directly 
converted  to  principal  stresses  and  directions  provided  that 
the  elastic  constants  for  the  material  are  known.  It  should 
be  noted  that  compressive  stress  gives  rise  to  extensional 
strain  relief  after  overcoring. 

The  situation  is  not  so  straight  forward  in  the  under- 
coring  case.  Here,  resort  must  be  made  to  elasticity  solutions 
that  express  the  stress  field  existing  before  trepanning  in 
terms  of  the  strain  relieved. 

Both  overcoring  and  undercoring  give  satisfactory 
results  when  used  in  macrostress  situations  provided  stress 
gradients  are  not  too  high,  and  are  the  basis  of  many  in 
situ  stress  measurement  devices.  Unfortunately,  there  are 
several  difficulties  associated  with  the  use  of  strain  relief 
methods  to  measure  residual  strains  in  relatively  small  samples 
^ock  or  total  strains  in  rock  masses  where  the  residual  part 
of  the  strain  relieved  is  a  significant  part  of  the  total. 

Successive  overcoring  of  the  same  strain  gages  shows 
that  strains  of  the  same  order  are  relieved  on  both  the  first 
and  the  second  cuts.  Experimental  results  are  reported  in  a 
later  section  for  a  medium  grained  granite.  If  one  or  more  of 
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the  microresidual  stress  mechanisms  are  acting,  all  of  the 
residual  strain  may  not  be  relieved  until  the  rock  is 
disaggregated  down  to  the  grain  size.  Even  then,  dislocations 
may  exist  in  the  crystal  lattice. 

If  the  strain  gage  size  is  of  the  same  order  as  the 
average  grain  size,  different  results  may  be  expected  from 
different  gages  with  the  same  orientation  depending  on 
whether  they  are  predominately  on  a  grain  showing  tensile 
or  compressive  relief.  On  the  other  hand,  if  the  gage  size 
is  many  times  the  average  grain  size,  more  consistent  results 
should  be  obtained. 

The  proximity  of  the  gages  to  the  relief  surface  may 
influence  the  results.  This  factor  is  generally  accounted 
for  in  trepanning  but  not  in  overcoring.  In  a  macro  field 
with  relatively  low  strain  gradients,  there  should  be  no 
problem.  Again  the  difficulty  is  associated  with  measuring 
microstrain  fields.  Even  in  macro  fields,  difficulty  can 
arise  if  a  plastic  zone  develops  due  to  the  stress  concentration 
around  the  hole . 

If  the  residual  strain  field  is  ideally  elastic,  there 
should  be  no  difficulty  with  time  dependent  effects.  However, 
time  dependent  effects  have  been  observed  during  some 
preliminary  tests  run  at  the  South  Dakota  School  of  Mines  and 
Technology.  Time  dependent  effects  have  also  been  reported 
by  Varnes  and  Emery  (loc.  cit.) 

Temperature  effects  previously  discussed  may  also 
effect  the  magnitude  of  residual  strains  relieved.  Further 
work  on  temperature  effects  is  planned. 

There  are  difficulties  associated  with  converting  strains 
relieved  into  stresses.  It  is  likely  that  the  material 
properties  at  this  scale  are  anisotropic  and  vary  from  grain 
to  grain.  Consequently,  the  use  of  large  scale  moduli  or  an 


isotropic  modulus  may  not  be  valid. 

Another  potential  difficulty  with  strain  relief 
measurements  occurs  if  a  strain  gage  spans  a  microcrack 
in  the  rock.  If  the  microcrack  opens  during  the  relief 
process,  seme  of  the  strain  recorded  by  the  gage  may 
actually  be  rigid  body  movement  and  not  deformation  of  the 
material 

Further  research  is  required  to  more  definitely 
establish  the  effects  of  the  above  difficulties.  Some 
recently  completed  work  is  reported  in  following  sections. 
Laboratory  Results 

Laboratory  tests  have  been  completed  on  a  section  of  6 
inch  core  of  Redfield  granite,  .1  5/8  inches  thick.  This 
granite  has  been  used  .in  previous  studies  of  residual  strain 
in  rock  and  has  been  described  as 

The  granite  ( Precambr ian )  from  Redfield, 

South  Dakota,  is  a  light  reddish  brown,  medium 
to  coarse-grained  rock  which  consists  of 
major  amounts  of  pink  orthoclase  feldspar,  grey 
vitreous  quartz  and  minor  amounts  of  black 
flakes  of  biotite.  In  thin  section,  the  rock 
consists  of  equigranular ,  anhedral  to  subbed ral 
minerals  of  fine  to  medium-grained  (0.3-3. 6  mm) 
quartz  (37%);  medium  to  coarse-grained  microcline 
(51%);  coarse-grained  (3.6mm)  plagioclnse  (10%); 
fine-grained  (0.5  mm)  biotite  (2%)  and  traces  of 
magnetite  and  apatite.  The  minerals  do  not  appear 
to  be  preferentially  oriented;  the  minerals  contain 
numerous  small  fractures  both  at  and  across  grain 
boundaries . 

The  section  of  core  used  had  a  smooth  surface  on  all  sides  . 


The  sample  size  and  strain  gage  locations  are  shown  on  Fig.  6 
The  strain  gages  were  initially  zeroed  and  allowed  to 
stabilize  under  the  cooling  water.  The  first  cut  was  made 
from  the  back  side  using  a  3  1/4  x  3  3/16  inch  bit.  This 
operation  overcored  the  inner  six  strain  gages  while 
simultaneously  undercoring  the  outer  five  strain  gages. 

Typical  data  is  shown  on  Fig.  7  for  the  rosette  formed 
by  the  three  1/4  inch  single  gages,  9,10,  and  11.  This 
rosette  was  overcored  by  drilling  from  the  back  side.  The 
innermost  rosette  formed  by  1/8  inch  gages  6, 7, and  8  was  also 
overcored  by  the  3  1/4  inch  bit.  As  shown  on  Fig.  7,  readings 
were  continued  after  the  bit  broke  through  in  order  to  check 
for  possible  time  dependent  effects.  Readings  were  discontinued 
when  two  successive  readings  showed  no  further  significient 
change  in  strain. 

The  initial  change  in  strain,  which  occured  when  a 
3/16  inch  deep  kerf  was  cut  on  the  opposite  side  of  the 
sample,  was  greater  for  the  inner  rosette.  This  is  contrary 
to  what  might  be  expected.  In  a  macrostress  field,  with 
small  strain  gradients,  one  would  expect  that  both  rosetts 
should  indicate  the  same  magnitude  of  strain  relief.  This 
is  taken  to  be  evidence  that  the  residual  field  in  this 
sample  is  micro  rather  than  macro  in  nature. 

It  is  interesting  to  note  that  the  principal  directions 
indicated  by  the  two  rosetts  are  relatively  consistent.  For 
the  inner  rosette,  the  principal  direction  was  found  to  be  135 
degrees  clockwise  from  gage  6  while  the  corresponding  direction 
from  the  next  outer  rosette  was  143  degrees  clockwise  from  gage 
9.  Consistency  in  directions  when  measuring  residual  strains  by 
relief  methods  has  previously  been  reported, 

Because  of  the  relatively  small  amount  of  strain  relieved 


Fig.  7:  Strain  relieved  by  cutting 
from  the  back  with  a  3  1/4"  bit. 
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at  kerf  depths  of  11/16,  15/16,  and  1  3/16  inches,  these 
readings  were  averaged  to  obtain  new  initial  strain  readings 
to  compare  with  strain  readings  obtained  after  the  bit  broke 
through  from  the  back  side.  In  this  case,  there  appears  to 
be  no  correlation  between  either  the  magnitude  of  the  strain 
releived  or  in  the  principal  directions  computed  from  these 
strains.  A  comparsion  of  the  results  from  the  two  rosettes  that 
were  overcored  is  shown  in  Table  1. 

After  the  above  test  had  been  completed,  the  same  sample 
was  cored  again.  This  time  with  a  15/16  x  7/8  inch  bit  from 
the  gage  side.  Thus  the  innermost  rosette  was  again  overcored 
while  the  next  outer  rosette  was  undercored.  The  strain 
changes  registered  by  gages  6,7  and  8  respectively  were-22, 

-49,  and  +  6.  The  principal  direction  was  found  to  be  46  degrees 
clockwise  from  gage  6.  Neither  the  magnitudes  nor  the  principal 
direction  ay  pear  to  correlate  with  the  overcoring  data  shown  in 
Table  1. 

In  both  overcoring  tests,  there  were  examples  of  strain 
gages  facing  each  other  across  the  kerf  of  the  drill  bit,  e.g., 
gages  6  and  9  in  the  second  test.  One  pair  of  these  gages 
in  each  test  registered  changes  in  strain  with  opposite  signs 
indicating  extension  on  one  side  and  contraction  on  the  other 
side  of  the  kerf.  This  is  taken  as  further  evidence  of  the  micro 
nature  of  the  residual  strain  field  in  this  piece  of  granite. 

A  macroresidual  strain  field  would  have  registered  consistent 
extension  or  contraction  on  both  sides  of  the  kerf  provided 
that  strain  gradients  were  relatively  small. 

further  confirm  the  micro  nature  of  the  residual 
strain  field  and  the  influence  of  gage  size  relative  to 
grain  size,  a  temperature  dependent  test  was  devised.  This 
test  was  conducted  on  the  7/8  inch  core  and  the  annular  shell 
remaining  from  the  second  overcoring  and  still  having  strain 


Table  1:  Comparison  of  results  from  simultaneously  over 
corinq  two  strain  gage  rosettes. 
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gages  6  through  11  intact.  The  inner  rosette  and  two  of 
the  gages  on  the  outer  shell  were  treated  as  live  gages 
while  the  remaining  gage  on  the  outer  shell  served  as  the 
temperature  compensating  gage.  Initial  readings  were  taken 
until  they  were  stable  at  room  temperature.  The  two  samples 
of  granite  were  then  placed  in  a  freezer  at  zero  degrees 
Fahrenheit.  Fig.  8  shows  the  change  in  strain  on  the  three 
active  gages  making  up  the  inner  rosette  when  gage  10  was 
used  as  the  common  temperature  compensation  gage.  The  same 
test  was  run  with  gage  number  11  on  the  outer  ring  acting 
as  the  compensating  gage. 

As  shown  on  Fig.  8,  temperature  compensation  was  not 
effective  throughout  the  duration  of  the  test.  The  immediate 
drop  in  strain  that  occured  when  the  samples  were  placed  in 
the  freezer  was  probably  due  to  surface  cooling  and  the  fact 
that  the  amount  of  surface  area  per  unit  volume  was  different 
for  the  7/8  inch  core  and  the  outer  shell  containing  the 
compensating  gage.  However,  if  temperature  compensation  were 
effective,  after  a  period  of  time,  the  gages  should  have 
again  registered  no  strain.  The  strain  change— time  curve 
shown  on  Fig.  8  indicates  that  after  80  minutes  at  zero 
degrees,  the  gages  registered  significant  strains  showing 
extension  at  gages  6  and  8  and  contraction  at  gage  7. 

Results  from  the  other  temperature  test  when  gage  11  was 
used  as  the  compensating  gage  were  similar  to  those  shown 
on  Fig.  8  except  that  the  magnitudes  of  the  strains  relieved 
were  315,  241,  and  313  ppm  respectively  for  gages  6,7  and  8 

after  approximately  95  minutes  at  zero  degrees.  The  principal 
directions  computed  from  the  strain  changes  were  26  degrees 
for  the  case  shown  on  Fig.  8  and  30  degrees  for  the  other 
case.  Again  the  magnitudes  of  the  strain  changes  do  not 
correlate  but  their  principal  directions  agree  reasonably 


well . 
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.  8:  Strain  change  verses  time  resulting 
from  temperature  change. 


These  results  may  be  interpreted  as  follows:  the  large 
difference  in  strain  marnitudes  between  the  two  tests  is 
probably  due  to  the  fact  that  the  strain  gages  used  for 
temperature  compensation  did  not  span  a  representative 
distribution  of  the  minerals  making  up  the  rock,  i.e.,  the 
gage  length  was  too  small  in  comparsion  to  the  mean  grain 
size.  In  one  case  the  grains  under  the  compensating  gage  had 
a  markedly  different  effective  coefficient  of  thermal 
expansion  than  in  the  other.  If  the  only  difference  between 
the  two  tests  wore  the  coefficient  of  thermal  expansion  of 
the  material  under  the  compensating  gage,  all  of  the  live 
gages  should  have  registered  the  same  change  in  strain. 

Since  this  was  not  the  case,  it  seems  likely  that  the 
residual  strain  field  is  temperature  dependent  in  this 
sample.  This  position  is  supported  by  the  relatively  good 
.agreement  in  the  principal  directions  computed  from  the 
test  data. 

Qualitative  System  Models 

The  models  shown  of  residual  stress  mechanisms  are 
readily  understood  when  they  are  considered  one  at  a  time 
as  previously  discussed.  However,  in  the  actua^.  case  it 
appears  likely  that  two  or  more  of  the  mechanisms  act 
simultaneously.  Furthermore,  a  large  number  of  self- 
equilibrating  mechanisms  interact  with  one  another.  The 
behavior  of  such  complex  systems  is  rather  difficult  to 
visualize.  As  a  consequence,  Varnes  has  fabricated  and 
tested  physical  models  to  aid  in  under stand ing  residual 
stress  phenomena.  These  models  consist  of  elastic, 
plastic,  and  viscous  elements  and  are  described  in 
Varnes  ( 19G9 ) . 
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As  an  alternative  to  physical  models,  numerical 
models  can  be  used.  Numerical  models  have  the  advantage 
of  being  easily  adapted  to  new  configurations  as  they  are 
encountered.  Also,  parameter  studies  can  be  made  to 
determine  the  sensitivity  of  the  behavior  to  changes  in 
relative  stiffness,  anisotropy,  viscosity,  etc. 

A  residual  stress  modeling  digital  computer  program 
is  presently  being  developed  at  the  South  Dakota  School 
of  Mines  and  Technology.  The  first  step  has  been  to 
develop  a  finite  element  program  that  employs  self-equilib¬ 
rating  elements.  An  example  of  such  an  element  is  shown 
on  Fig.  9.  The  element  is  composed  of  eight  separate 
elastic  springs  connected  at  five  nodal  points.  The 
interior  nodal  point  may  be  eliminated  by  a  substructure 
analysis  leaving  only  the  four  corner  nodes.  Each  corner 
node  has  two  degrees  of  freedom.  This  basic  square  or 
rectangular  rock  sample.  An  example  is  shown  on  Fig.  9. 

Each  element  in  the  rectangular  array  of  elements 
(the  model)  shown  on  Fig.  9  is  connected  to  its  adjoining 
elements  only  at  the  corner  nodal  points.  Residual 
stresses  are  induced  into  the  elements  by  allowing  the 
inner  springs  in  each  element  to  have  a  different  coefficient 
of  thermal  expansion  than  the  springs  forming  the  outer 
square.  The  temperature  of  the  model  is  then  changed 
uniformly.  If  the  inner  springs  have  a  higher  coefficient 
of  thermal  expansion  and  if  the  temperature  change  is  an 
increase,  the  inner  springs  in  each  element  will  develop 
compression  and  the  outer  springs  tension.  The  model  at. 
this  stage  contains  residual  stresses. 

A  cut  in  an  outer  spring  of  an  element  in  the  model 
may  be  simulated  by  reducing  the  stiffness  of  that  partic¬ 
ular  spring.  It  should  be  noted  that  reducing  the  stiffness 
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of  a  particular  spring  in  one  of  the  outer  elements 
results  in  a  change  in  displacement  of  all  of  the  nodes 
in  the  model. 

Some  preliminary  results  obtained  from  the  finite 
element  model  are  shown  on  Fig  .10  where  the  change  in 
length  of  two  typical  springs  from  the  model  are  plotted 
versus  the  percent  reduction  in  stiffness  of  the  outer 
horizontal  spring  connecting  nodes  3  and  4  in  element 
number  3.  The  parameter  "c"  equals  the  ratio  of  the 
stiffness  of  the  inner  springs  in  an  element  to  the 
stiffness  of  the  outer  springs  in  the  element.  All 
elements  are  identical  prior  to  the  reduction  in  stiffness 
of  the  spring  connecting  nodes  3  and  4.  The  interior 
coefficients  of  thermal  expansion  are  1.414  times  the 
outer  coefficients. 

The  results  shown  on  Fig. 10  illustrate  the  importance 
•of  the  stiffness  ratio  c  to  the  response  of  the  system  to 
the  reduction  in  stiffness  of  one  outer  spring.  Cutting 
the  outer  spring  between  3  and  4  results  in  an  extentional 
change  in  length  on  both  the  top  and  bottom  surfaces  of 
the  model  when  c  -  4  and  a  contract ional  change  when 
c  =  0.25  up  to  the  point  where  the  cut  has  reduced  the 
stiffness  by  75%.  At  this  point,  the  spring  connecting 
nodes  4  and  5  had  been  dominated  by  lateral  contraction, 
after  this  point,  vertical  movement  causing  extension 
begins  to  play  an  important  role  in  the  change  in  length 
of  the  spring.  The  single  most  impressive  feature  of  the 
results  shown  on  Fig. 10  is  the  complexity  in  the  response 
of  such  a  simple  model,  it  should  be  noted  that  the  change 
m  strain  resulting  from  the  first  cut  on  the  back  side 
of  the  granite  sai  pie  discussed  in  the  last  section 
might  be  explained  by  a  numerica]  model  similar  to  the 
one  discussed  here. 


Change  in  strain 
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Fig.  .10:  Results  from  numerical  model. 


It  is  easy  to  imagine  more  complex  models.  For 
example,  the  simple  square  element  could  be  replac'd 
by  a  self-equilibrating  hexagon  or  octagon  with  interior 
springs  radiating  from  a  single  interior  nodal  point  and 
each  having  its  own  chacteristics .  Friction  elements  and/ 
or  vicous  elements  could  be  added  in  series  and/or  parallel 
to  the  spring  elements.  A  whole  realm  of  mixed  element 
models  could  be  developed  having  different  shapes  and 
chacteristics  to  attempt  to  better  simulate  the  behavior 
of  actual  rocks.  All  of  the  suggested  models  could  be 
extended  to  three  dimensions.  T^r-tically  speaking,  storage 
space  in  the  computer  soon  becomes  a  problem  and  limits 
extent  of  the  simulation. 

At  this  stage  of  development,  the  primary  value  of 
both  physical  and  numerical  models  appears  to  be  their 
use  as  an  aid  in  understanding  the  complex  behavior  of 
large  systems  containing  self-equilibrating  elements. 
Presently,  it  is  not  feasible  to  obtain  quantitative 
information  from  these  models,  hence  they  are  referred  to 
as  qualitative  models. 

Discussion 

Results  from  both  laboratory  tests  and  a  very  simple 
mathematical  model  indicate  that  the  response  of  arrays 
of  interconnected  self-equilibrating  units  is  complex  and 
not  easily  visualized.  The  laboratory  test  results  seem 
to  indicate  that  at  least  part  of  the  residual  strain 
field  is  temperature  dependent  and  the  size  of  the  strain 
gage  relative  to  the  grain  size  is  probably  an  important 
parameter . 

The  likely  existence  of  a  temperature  dependent 
residual  strain  field  suggests  that  thermal  weakening  of 
some  rocks  may  be  due  to  residual  stresses  changing  with 
temperature.  It  is  anticipated  that  a  set  of  controlled 
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experiments  will  be  devised  to  investigate  this  hypothesis. 

The  fact  that  both  the  laboratory  results  and  the 
numerical  model  indicate  that  strains  may  be  registered  on 
the  opposite  side  of  a  sample  from  where  a  cut  is  made. 

It  appears  likely  that  there  may  be  a  shape  effect,  i.e., 
the  size  and  shape  of  the  sample  probably  influence  the 
amount  and  distribution  of  the  strain  relieved. 

Future  work  on  residual  stresses  in  rock  will  entail 
studies  of  the  gage  size  to  grain  size  ratio  as  well  as 
the  grain  size  to  sample  size  ratio.  More  tests  will  be  run 
to  attempt  to  define  temperature  dependence  of  the  residual 
strain  fields  in  rock.  Strains  relieved  at  different 
positions  on  the  surface  of  a  sample  due  to  a  hole  cut  at 
one  point  should  be  measured.  Numerical  modeling  will 
continue  to  help  in  the  understanding  of  residual  stresses 
in  rock. 

V  Data  Reduction 

Reduction  of  measured  strain  data  from  the  deep  hole 
device  to  in  situ  stresses  can  be  accomplished  in  any  of 
several  ways.  A  combination  analytical  -  experimental 
technique  was  described  in  the  final  report  of  the  first 
year's  work  on  this  project  and  also  in  a  paper  by  Hoskins 
and  Oshier  published  in  the  14th  Symposium  on  Rock  Mechanics. 
This  method  is  summarized  in  the  section  entitled  Experimen¬ 
tal  Method. 

Sateesha  working  on  this  project  developed  a  three- 
dimensional  finite-element  code  for  reduction  of  data  from 
the  deep  hole  device.  Portions  of  his  thesis  are  summarized 
in  the  section  entitled  Numerical  Method. 

Russell  also  working  on  this  project  produced  an 
approximate  analytical  solution  for  data  reduction.  His 
work  was  presented  at  the  15th  Symposium  on  Rock  Mechanics 
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and  is  summarized  in  the  section  entitled  Analytical  Method. 
Experimental  Method 

Hiramatsu  and  Oka  (1962)  have  given  the  components  of 
stress  in  an  elastic  infinite  body  surrounding  a  cylindrical 
borehole  as  follows: 
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In  the  sidewalls  of  a  borehole  subjected  to  a  fluid 
pressure  P 
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When  r,  9  and  5  are  the  axes  of  a  cylindrical  coordinate 
system,  IT*  =  (Poisson's  ratio)  and  the  ot  ’  S 

and  Y's  are  combinations  of  direction  cosines. 
The  stress  strain  relations  can  be  written 

^E=l,rV'(  %  **i  > 

£$  E  =  -  V(  cfr  +  crQ) 

q  G  ~  Q 

Now  from  the  three  strain  gage  rosettes  in  the  device 
we  shall  determine  f  ^  ,  and  VV  Q  each  at  (ft ,  <f>4->20 

and  (P  +  24-0  in  the  borehole  walls.  We  thus  have  9 
equations  with  which  to  calculate  the  3  principal  stresses 
Pi.  P2<  and  P3  and  their  directions  of  action  relative 
to  the  borehole  axis  and  an  arbitrarily  selected  direction 
This  much  of  the  analysis  is  similar  to  Leeman  and  Hayes 
(1966) .  The  trepanning  technique  does  not  accomplish 
complete  stress  relief,  however.  The  percentage  of  stress 
relief  achieved  depends  upon  the  dimensions  of  the  strain 
gages  and  the  trepanning  hole  and  theposition  of  the  gages 
relative  to  the  hole.  In  addition  the  friction  bonded 
strain  gages  may  not  be  100%  efficient  at  responding  to 
changes  in  strain.  These  two  factors  can  be  combined  and 
experimentally  determined  to  yield  an  overall  value  for 
the  efficiency  of  this  measurement  device. 

The  efficiency  of  friction  bonded  strain  gages  them¬ 
selves  is  a  function  of  their  design  and  the  normal  stress 
used  to  hold  them  in  contact  with  the  borehole  walls. 
Laboratory  tests  so  far  show  that  85%  is  a  typical  figure 
for  the  efficiency  of  the  gages  that  we  have  made.  We 
plan  initially  at  least  to  individually  calibrate  each 


O  e  e  o  friction  strain  gages  before  using  them  in  the 
deep  hole  device. 

Numerical  Method 

General 

The  Finite-Element  Method,  now  widely  accepted  in 
engineering  mechanics,  is  particularly  well  suited  for  the 
solution  of  boundary  value  problems  of  a  geological  nature 
It  can  be  easily  adapted  to  the  solution  of  systems 
characterized  by  inelastic,  anisotropic,  heterogeneous 
material  properties,  of  any  exterior  configurations 
containing  structural  discontinuities,  and  subjected  to 
any  statically  viable  boundary  conditions  of  load  or 
displacement.  The  main  advantages  of  the  finite-element 
method,  as  compared  to  other  numerical  techniques,  such  as 
the  classical  finite  difference  formulation, have  been 
described  by  Zienkiewicz  (1971),  Clough  (1965),  and  Felippa 
and  Clough  (1970).  Exclusive  emphasis  in  this  chapter  will 
therefore  concern  the  applications  of  finite-element 
methods  for  solving  geological  problems,  and  an  important 
practical  problem,  for  which  a  closed  form  solution  is 
available,  will  be  presented  to  illustrate  the  method. 

The  range  of  specific  problems  encountered  in  geology 
is  large,  but  some  typical  examples  that  have  been  fully 
worked  out  utilizing  the  finite-element  method  include 

(I)  behavior  of  rocks  under  static  and  dynamic  loads, 

(II)  seismic  response  and  propagation  of  stress  waves, 

(HI)  heat  conduction,  (IV)  fluid  flow  in  porous  media, 

and  (V)  distribution  of  magnetic  and  graviational  potential, 
Many  of  the  fundamental  problems  in  structural  geology, 
geophysics,  geohydrology,  geomorphology,  glacialogy,  and 
engineering  geology  can  be  included  in  the  above  categories. 

Very  few  authors  have  attempted  to  solve  geological 
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problems  applying  the  principles  of  classical  mathematical 
physics.  Hafner  (1951)  used  the  Airy  stress  function 
of  classical  elasticity  to  obtain  two-dimensional  solutions 
for  stress  distributions  caused  by  several  arbitrarily 
defined  forms  of  boundary  forces,  and  constructed  the 
potential  fault  surfaces  based  on  the  original  stress 
distribution.  Sanford  (1959)  used  an  approach  similar 
to  that  of  Hafner  for  the  theoretical  analysis  of  fault 
structures.  A  Fourier  series  approximation  was  used  to 
define  the  displacement  boundary  condition  at  the  base 
of  a  homogeneous  elastic  layer.  Howard  (1966)  applied 
Sanford's  theory  for  the  analysis  of  William  Range  Thrust 
Fault  in  Middle  Park,  Colorado. 

The  stress  function  solutions  thus  obtained  from  highly 
idealized  models  have  proven  to  be  both  stimulating  and 
valuable.  However,  with  few  exceptions  modification  of 
stress  function  solutions  to  better  fit  given  field 
situations  has  rarely  been  attempted.  Even  with  the  avail¬ 
able  solutions,  the  restrictive  assumptions  involved  are 
often  ignored.  To  start  with,  our  conception  of  the  real 
problem  is  in  itself  a  model.  In  geologic  problems 
isotropic  elasticity  cannot  always  be  assumed,  body  forces 
must  be  taken  into  account  and  the  lateral  and  vertical 
variation  in  material  properties  that  is  generally  present 
must  be  considered.  Well-defined  discontinuities  may  be 
present  and  progressive  failure  of  some  kind  may  be  induced 
by  the  assumed  boundary  force  and  displacement  conditions. 
Tie  situation  gets  more  complicated  if  the  mechanism  is 
path  dependent,  i.e.,  if  the  final  state  of  stress, 
strain  or  displacement  depends  upon  the  manner  in  which 
the  surface  tractions,  displacements  and  body  forces 
obtain  their  final  values. 


44 


The  effect  of  the  complexities  cited  above  cannot  be 
adequately  determined  without  recourse  to  analytical  and 
experimental  evaluation.  Existing  closed-form  mathematical 
tools  may  be  suitable  for  solution  of  only  a  few  such 
problems,  depending  on  the  nature  of  the  parameter  (e.g. 
some  classes  of  anisotropy).  In  each  of  the  above  situations, 
the  finite-element  method  is  sufficiently  versatile  for 
definite  evaluation  of  the  problem. 

To  geologists,  the  finite-element  method  has  thus 
opened  up  the  possible  exploration  of  stress  distributions 
in  a  large  variety  of  complex  situations  of  geologic  impor¬ 
tance.  Such  distributions  are  important  not  only  for 
problems  of  faulting  but  also  for  folding  and  ductile 
deformation  in  general.  Dietrich  and  Carter  (1968)  have 
used  the  finite-element  method  to  study  the  stress  history 
of  folding . 

The  deformation  of  rock  associated  with  cylindrical 
underground  excavation  has  been  investigated  by  several 
authors.  With  few  exceptions  the  problem  has  been  treated 
as  a  two-dimensional  case  and  one  of  the  principal  stress 
directions  is  frequently  assumed  to  coincide  with  the  axis 
of  the  borehole.  The  orientation  of  the  principal  stresses 
underground  is,  in  general,  not  known  and  it  must  be 
assumed  that  the  axis  of  a  borehole  does  not  coincide  with 
a  principal  stress  direction.  Such  a  general  situation  can 
only  be  analyzed  by  a  three-dimensional  analysis  and  this 
provides  an  excellent  testing  ground  to  check  the  validity 
of  three-dimensional  finite-element  analysis. 
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THREE-DIMENSIONAL  FINITE-ELEMENT  CTRESS  ANALYSIS  COMPUTEP 
MODEL  STUDIES 

The  finite-element  technique  used  in  this  study 
is  based  on  the  displacement  or  stiffness  method  of 
analysis  first  developed  by  structural  engineers  for  air¬ 
craft  industry.  These  programs  are  capable  of 
analyzing  any  three-dimensional  structure  subjected 
to  any  general  type  of  loading.  However,  the 
applied  stress  field  is  unknown  in  the  borehole¬ 
trepanning  hole  intersection  problem  under  consideration 
and  paradoxically,  knowing  the  stress  field  acting 
on  the  model  is  tantamount  to  determining  the  primary 
stress  field.  To  circumvent  this  difficulty,  the 
method  of  strain  coefficients  is  adopted  and  is 
described  in  the  following  paragraphs. 

The  method  of  strain  coefficients  has  been  used 
by  Clough  to  solve  a  two-dimensional  tunnel  problem 
where  the  cross  section  was  not  circular.  For  the 
borehole-trepanning  hole  intersection  problem,  the  following 
can  be  considered  a  general  three-dimensional 
algorithm  for  determining  the  primary  stresses  from 
strain  measurements  on  the  sidewalls  of  a  borehole. 


1. 


Choose  a  three-dimensional  finite-element  model 


consisting  of  a  vertical  borehole  in  the  middle  as  shown  in 
Figure  11  The  model  consists  of  a  finite  number  of 
elements  interconnected  at  corners  or  nodal  points.  This 
physical  idealization  of  a  solid  as  an  assemblage  of  finite- 
elements  involves  no  mathematical  approximation.  The  only 
approximation  is  in  assuming  a  displacement  function  for 
each  element.  The  usual  procedure  is  to  use  a  linear  dis¬ 
placement  function.  The  finite-element  solution  can  be 
shown  to  converge  to  the  exact  solution  as  element  size  is 
reduced  to  a  point. 

2.  Apply  six  unit  stress  fields  (a  o  a  x  x 

x  y  z  xy  yz 

Tzx^  separately  on  the  model  without  trepanning  holes  and 
determine  the  strains  at  gage  points  from  a  three- 


dimensional 

finite' 

-element 

analysis . 

Let 

these  strains 

(Figure  12 

be 
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G  =  11tt/6 

(fig.  3.1) 

. ,6,  where  the  first  subscript  ’ b ’  indicates 

strain  in  the  model  before  trepanning  and  the  second  sub¬ 
script  gives  the  gage  location.  The  superscript  (i) 
represents  the  applied  stress  field.  Thus  i=l  corresponds 
to  i=2  corresponds  to  etc. 

3.  Now  consider  the  model  with  three  trepanning  holes 
120°  apart  on  the  sidewalls  of  the  borehole.  Again  apply 


Applied  unit  stresses  and  associated  strains 
at  gage  points  for  the  model  without  holes 


six  unit  stress  fields  (o  a  a  t  t  t  ) 

x  y  z  xy  yz  zx 

separately  and  determine  the  strains  at  the  gage  points  from 
a  three-dimensional  finite-element  analysis.  Let  these 
strains  (Figure  13)  be 
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i=l/2, . ,6,  where  the  first  subscript  ' t 1  indicates 

strain  in  the  model  with  trepanning  holes. 

4.  Compute  the  resulting  relaxation  strains  or  strain 
coefficients  by  subtracting  strains  obtained  in  step  3  from 
.  the  corresponding  strains  determined  in  stop  2.  Thus 


-j  =  i  2  9-  i=]  2 

rj  bj  ctj  '  3  1,Zr . 1  L,Z 


6, 


where  the  subscript  ’r*  indicates  the  relaxation  strain. 

5.  Let  the  measured  value  of  relaxation  strains  in  a 
three-dimensional  stress  field  be  denoted  by 
e 
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Cm2 

e  _  at 
m3 

position 
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-  0=ir/2 , 

m4 
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e  r  at 
m6 

position 
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-  0=771/6, 

m7 
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(3) 

-  0  =  1171/6, 

and  e 

where  the  first  subscript  'm'  indicates  the  Pleasured  strains 
and  the  second  subscript  gives  the  gage  location. 

Let  us  assume  that  the  strain  relieved  at  any  of  the 


three  6  positions  does  not  affect  the  measured  strains  at 
the  remaining  two  positions.  This  assumption  is  valid 
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since  the  trepanning  hole  diameter  is  very  small  in 
comparison  with  borehole  diameter  and  the  spacing  between 
any  two  holes  is  several  times  the  hole  diameter  along  the 
borehole  circumference.  Let  us  further  assume  that  the 
model  material  behaves  in  a  linearly  elastic  manner  and 
the  principle  of  superposition  is  valid.  The  relationship 
between  the  measured  strain  and  the  unknown  stress  com¬ 
ponents  can  now  be  written  by  summing  the  products  of  the 
stress  components  and  the  associated  strain  coefficients 
for  each  gage  location.  Therefore 

e(1)c  +  c(2)o  +  e(3)o  +  e(4)x  +  c(5>t  +  s(6>t 

ri  x  rj  y  rj  z  r;j  xy  rj  yz  rj  zx 


emj'  j  lf2' 


or  (e  ) to}  =  (em) 


(3.24) 

(3.24a) 


Since  there  are  nine  equations  and  only  six  unknowns, 
the  system  is  overdetermined  and  the  additional  measure¬ 
ments  are  said  to  be  overcomplete.  Such  cases  are  very 
common  in  practice  since  we  measure  more  quantities  than 
there  are  unknowns,  to  avoid  the  influence  of  statistical 
fluctuations.  The  problem  can  be  solved  by  two  different 
approaches.  In  the  first,  six  independent  equations  are 
arbitrarily  chosen  and  solved  for  the  six  unknown  stress 
components.  All  possible  combinations  of  six  equations 
that  do  not  result  in  a  singular  system  are  then  solved. 
The  solution  vectors  so  obtained  are  then  used  for  com¬ 
parison  and  cross  checking.  Using  equations  (3.19)  and 
(3.20),  it  can  be  shown  that  only  18  combinations  are 
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are  possible.  Many  of  these  combinations  may  result  in  ill- 
j  .  . 

conditioned  system  of  equations.  In  the  second  approach, 
which  is  more  rational,  we  solve  the  problem  by  means  of  the 
maximum  likelihood  method  from  the  theory  of  adjustments. 

) 

The  best  estimates  for  the  unknown  stress  components  are 
obtained  by  applying  the  method  of  least  squares.  The 
normal  equations  are  obtained  as 

(er)T* (W) • (er) *{o}  =  (er)T- (W) *{om}  (3.24b) 

where  (W)  =  weight  matrix,  and  the  superscript  T  designates 
the  transpose  of  a  matrix. 

6.  Compute  the  principal  stresses  o^,  °(2)  and  °(3)' 

and  their  orientations  from  the  six  stress  components 
obtained  in  step  5. 

The  primary  stresses  in  the  rock  are  now  completely 
determined. 

Several  computer  programs  were  developed  to  carry  out 
the  analysis  described  in  the  above  steps.  These  programs 
are  completely  compatible  with  one  another  and  a  description 
of  their  use  is  given  in  the  Appendix. 

A  large  scale,  high-speed  digital  computer  is  required, 
owing  to  the  large  number  of  elements  normally  required  to 
adequately  model  the  borehole-trepanning  hole  intersection. 

A  CDC-3400  computer  with  32K  word  storage  was  used  in  this 
study. 

SUBSTRUCTURE  METHOD  OF  ANALYSIS 

The  number  of  finite-elements  required  to  obtain  the 
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stress  distribution  at  the  borehole- trepanning  hole  inter¬ 
section  is  so  large  that  the  available  computer  facilities 
cannot  conveniently  accomodate  the  problem.  To  overcome 
this  difficulty,  the  SUBSTRUCTURE  METHOD  of  analysis  is 
employed . 

The  substructure  method  of  analysis  is  well  known  to 
structural  engineers  and  is  usually  applied  to  such  large 
structures  as  aircraft  frames,  multi-story  buildings  and 
ships.  Each  of  these  can  be  considered  to  consist  of  a 
number  of  substructures  obtained  from  structural  part¬ 
itioning.  Even  though  the  partitioning  can  be  arbitrary, 
it  is  preferable  to  make  structural  partitioning  correspond 
to  physical  partitioning.  An  aircraft  frame,  for  example, 
might  be  considered  to  consist  of  wing,  tail  and  fuselage 
components . 

The  fundamental  principles  of  substructure  analysis 
have  been  clearly  outlined  by  several  authors,  for  example 
Przemien ie^c^i  1968  using  a  displacement  approach  and 
Argyris  /  using  redundant  interaction  force  concepts.  The 
most  compact  treatment,  and  that  most  generally  suitable 
for  programming  is  the  displacement  method  and  this  is  used 
in  the  following  condensed  theory. 

In  the  substructure  method,  the  substructures  are 
treated  as  if  they  were  complicated  finite-elements,  inter¬ 
connected  at  terminal  nodes  to  form  the  overall  structure. 

In  the  problem  under  consideration,  the  finite-element 
containing  the  trepanning  hole  is  treated  as  a  substructure. 
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The  stiffness  matrix  of  the  substructure  is  determined  by 
subdividing  the  substructure  into  a  number  of  smaller, 
simple,  finite-elements,  computing  the  simple  element  stiff¬ 
nesses  and  assembling  the  global  stiffness  matrix  of  the 
substructure,  and  condensing  the  substructure  stiffness  to 
eliminate  the  internal  degrees  of  freedom.  The  resulting 
substructure  stiffness  is  used  to  develop  the  overall 
stiffness  matrix  of  the  present  structure  which  is  analyzed 
for  any  applied  external  loading.  The  detailed  solutions 
for  the  substructure  are  obtained  by  determining  the 
internal  node  displacements  and  then  substituting  the  dis¬ 
placement  values  in  the  expressions  for  substructure 
stresses  and  strains. 

In  carrying  out  the  above  analysis,  the  continuity  of 
displacements  is  violated  at  the  substructure-parent 
structure  boundaries.  The  continuity  is,  of  course, 
satisfied  at  the  corner  nodes  of  parent  structure  element 
containing  the  substructure  by  virtue  of  the  application 
of  assembly  rules.  Since  the  interface  compatibility  is 
not  satisfied,  a  rigorous  mathematical  proof  of  an  upper 
bound  to  the  stiffness  cannot  be  obtained.  Nevertheless, 
the  above  procedure  is  being  used  increasingly  because 
it  usually  provides  useful  engineering  solutions. 

NUMERICAL  EXAMPLE 

Consider  a  three-dimensional  finite-element  model 
of  dimensions  24"x24"x36"  as  shown  in  Figure  14A .  Because  of 
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symmetry,  it  is  only  necessary  to  consider  an  octant  of  the 
model  as  shown  in  Figure  14b.  The  geometry  and  material 
properties  are  given  as  follows: 

Borehole  diameter  =  6  5/8" 

Trepanning  hole  diameter  =  1/4" 

Young's  modulus  =  8.57x  106  psi 

Poisson's  ratio  =  0.191 

The  above  geometry  and  material  properties  were  chosen 
because  they  represent  the  rock  block  models  of  Redwood 
granite  used  to  verify  the  proposed  method  of  analysis. 

The  parent  structure  (Figure  15,  is  modeled  into  an 
assembly  of  203  elements  and  the  substructure  (Figure  16) 
is  further  subdivided  into  49  elements  inter-connected  at 
130  nodal  points.  Nodes  3,  4,  5,  8,  9  and  10  are  common 
to  both  parent  and  substructure.  Six  types  of  loading  are 
applied  to  the  parent  structure.  The  strains  at  three  gage 
points  (Figure  17 )  around  the  trepanning  hole  are  determined 
from  substructure  finite-element  analysis  and  are  summarized 
in  Table  2 


Fig.  15  The  parent  structure 
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TABLE  2  STRAINS  AROUND  TREPANNING  HOLE 


EXTERNAL 

LOAD 

(psi) 

STRAINS  (f 
GAGE  1  t 

(Vertical) 

yin/in.)  AROUND 
GAGE  2 
(+120°) 

TREPANNING  HOLE 

GAGE  3 
(-120°) 

a  =+l 

X 

.01551 

.1345 

.1345 

0  =+l 
y 

-.00517 

-.0388 

-.0388 

a  =+l 
z 

.05208 

.01654 

.01654 

T  =  +  l 
xy 

.  2266 

-.3277 

.158) 

T  =  +  l 
yz 

-.011987 

.0497 

.0071 

T  =  +  l 
zx 

.0092 

.0156 

.2317 

The 

strains  at  these  three  gage  points  in  the  absence 

of  trepanning  holes  are 

obtained  from  a 

straight-forward. 

three-dimensional  finite 

-element  analysis  and  are  given  in 

Table  3 

TABLE 

3  STRAINS  ON  THE  SIDEWALLS 

OF  BOREHOLE 

EXTERNAL 

LOAD 

(psi) 

STRAINS  (-c 
GAGE  1  b 

(Vertical) 

uin/in.)  ON  THE 
GAGE  2 
(+120°) 

:  BOREHOLE  WALLS 

GAGE  3 
(-120°) 

a  =+l 

X 

-.07399 

.1345 

.1345 

0  =+l 

y 

.02287 

-.08181 

-.08181 

a  =+l 
z 

.1167 

.01202 

. 01202 

X  =  +  l 
xy 

.0231 

-.08483 

-.08483 

x  =+l 
yz 

0 

.03141 

-.03141 

x  =+l 
zx 

0 

-0.2407 

+0.2407 

The  relaxation  strains  are  now  obtained  as  the  differ¬ 
ences  between  the  values  listed  in  Tables  2  and  3  and 
are  presented  in  Table  .  4 >  The  three  observation  equations 
may  now  be  written  as 
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c<?> 

rl  rl 

r ( 6 )  (1) 
rl  x 

o 

Cm2 

f  (1)  A2)  (3) 
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e(*>  a 

r  2  y 

em3 

r(D  _(2>  (?) 
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t  (3.25) 

xy 

or  { e  } 

m 

(1)  =(er}(1){a} 

T 

yz 

T 

zx 

(3.26) 

TABLE 

4  RELAXATION  STRAINS  DUE  TO  TREPANNING 

EXTERNAL 

LOAD 

(psi) 

RELAXATION  STi  \INS  (c  uin/in 
GAGE  1  GAGfi  2 

(Vertical)  (+120°) 

.)  DUE  TO  TREPANG. 
GAGE  3 
(-120°) 

a  =1 

X 

.0895 

-.1376 

-.1376 

a  =1 

y 

-.028 

.043 

.043 

a  =1 
z 

-.06462 

.004516 

.004516 

T  =1 
xy 

.  2035 

-.2429 

-.2429 

o 

Tyz~  1 

-.011987 

. 01829 

.03851 

Tzx=  1 

.0092 

.2563 

-.0091 

“\ 

where  e  , , 
ml 

e  c  *>  are  the 

m2  m3 

relaxation  strains  measured  at 

IJ 

position  (1)  along  z  and  +120 

O  i  • 

directions . 

Two  more 

sets  of  equations 

of  the  above 

type  can  now  be 

obtained  for  positions  0=7tt/6  and  e=lln/6  by  a  simple 
transformation  of  components  of  stress.  For  the  systems  of 
axes  OX  Y'Z'  ano  0X"Y"Z"  rotated  through  120°  in  the  counter¬ 
clockwise  and  clockwise  directions  respectively  about 

Z-axis , 
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(3.27) 

(3.28) 


{ e  } <2)  =  (e  }  (2) {o} '  =  (e  }(2) (T) (2){o) 
m  r  r 

and  (e  }(3)  =  {e  }(3){o}"  =  (e  } (3) (T)  (3) (o } 
m  r  r 

where  {e  }  ^  and  {e  }  (^are  the  relaxation  strain  vectors 
m  m 

measured  at  8=7tt/6  and  0=11tt/6  positions  respectively,  and 
{a}'  and  {a}"  are  the  stress  components  in  OX'Y'Z1  and 
OX" Y" Z "  systems  respectively.  (T)  ^  and  (T)  ^  are  the 
coordinate  transformation  matrices  and  it  can  be  shown  that 
(10) 


(T) (2^ 


l2 

1 

m2 

« 

n2 

1 

21  jmi 

2m  1  n  1 

2n  1 1 1 

1 

l2 

1 

n2 

21 2m2 

2m  2  n  2 

2n  2 1 2 

2 

2 

2 

l2 

m2 

n2 

21  m 

2m  n 

2n  1 

3 

3 

3 

3  3 

3  3 

3  3 

1  1  1  2 

mim2 

n  in2 

(1 im2+l2mi ) 

(m  1  n2+in2n  1 ) 

(n ! l2+n2l , ) 

1  2  1  3 

in  2in  3 

n2n3 

(l2m3+l 3m2 ) 

(m2n  3+m  3n2 ) 

(n  2 1 3+n  3 12 ) 

l3ll 

ni3ini 

n3n  1 

( 1 3m 1 +1 im  3 ) 

(m3n  1  +m  1  r.  3 ) 

(n  3 1 1 +n 1 1 3 ) 

(3.29) 

where  the  direction  cosines  1,  m,  n  between  the  primed  and 
unprimed  coordinate  axes  are  defined  in  the  following  table: 


X 

Y 

Z 

X' 

li 

mj 

n  1 

Y' 

I2 

m  2 

n  2 

Z' 

I3 

m3 

n  3 

From  the  above  table  we  see  that  li=Cos  X'X,  m?=Cos  Y'Y  etc. 
An  equation  similar  to  (3.29)  may  be  written  for  (T)  ^ 3 ^ 
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For  the  problem  under  consideration,  the  transformation 
tables  of  direction  cosines  are 


X 

Y 

Z 

X' 

-1/2 

/3/2 

0 

Y* 

-  /3/2 

-1/2 

0 

Z' 

0 

0 

1 

and 


X 

Y 

Z 

X" 

-1/2 

-/3/2 

0 

Y" 

/3/2 

-1/2 

0 

Z" 

0 

0 

1 

Substituting  these  values  in  the  expressions  for 
and  (T) ^  and  combining  equations  (3.26),  (3.27)  and 
(3.28),  we  obtain 


(T) 


(2) 


{e  }  =  (A) {a}  (3.30) 

m 

where  {e  }  is  a  9x1  vector  corresponding  to  measured  relax- 
m 

ation  strains  at  three  0-positions  and  (A)  is  given  by 


.0895  - 

.028 

-.06462 

.2035 

-.1199 

.0092 

-.1376 

.043 

.004510 

-.  2429 

.0183 

.  2563 

-.1376 

.  043 

. 004516 

.  2429 

.0385 

-.0091 

.0895  - 

.028 

-.06462 

-.2035 

.0679 

.0992 

-.1073 

.0127 

.004516 

.2778 

.2128 

-.1441 

.1030  - 

.1976 

.004516 

.0349 

-.0271 

-.0288 

.0895  - 

.028 

-.06462 

.0016 

.052 

-.1084 

.1030  - 

.1976 

.004516 

-.0349 

.2311 

-.1123 

-.1073 

.0127 

.004516 

-.2778 

-.1146 

.0378 

Thus  nine 

equations  are 

obtained 

to  dote 

rmine  th 

a 
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unknown  stress  components.  From  this  overdetermined  set  of 
equations  one  can  solve  for  the  stresses  {a}  by  either  of 
the  two  methods  already  described  in  the  previous  sections. 
The  set  of  equations  (3.30)  thus  represent  the  necessary 
equations  to  determine  the  stress  tensor  in  any  linearly 
elastic ,  isotropic  and  homogeneous  rock  at  any  depth  in  the 
earth's  crust  from  the  strain  measurements  on  the  sidewalls 
of  deep  boreholes. 

Following  the  development  of  a  new  three-dimensional 
method  for  the  determine  U  a  of  primary  stresses  in  rock, 
an  experimental  investigation  has  been  planned  to  verify 
some  of  the  results  of  the  analysis.  The  experimental 

program  is  presented  in  the  next  section. 

O 

o 

o 
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EXPERIMENTAL  STUDIES 

INTRODUCTION 

A  program  was  established  for  testing  large  rock  block 
models  whereby  the  maximum  amount  of  information  could  be 
obtained  from  the  minimum  number  of  tests.  The  program  was 
developed  to  satisfy  two  specific  goals: 

(1)  To  verify  the  adequacy  of  the  model  components  that 
had  been  developed;  namely  model  material,  loading  apparatus 
and  instrumentation. 

(2)  To  yield  information  on  the  response  of  a  limited 
number  of  models  at  the  borehole-trepanning  hole  inter¬ 
section  so  that  the  results  may  be  compared  with  the 
numerical  solution  obtained  in  the  previous  chapter. 

Two  kinds  of  tests  were  carried  out  on  two  different 
rectangular  block  models  of  the  same  material.  The  differ¬ 
ence  between  the  two  models  results  from  different  borehole 
orientation  and  loading  with  reference  to  borehole  axis. 
Nevertheless  exactly  the  same  scheme  and  sequence  of  loading 
were  used  for  both  models. 

Strain  gage  rosettes  were  placed  at  selected  locations 
on  the  sidewall  of  the  borehole.  The  rosettes  were  mounted 
on  the  blocks  so  that  one  strain  gage  was  parallel  (or 
nearly  parallel)  to  the  borehole  axis.  This  selected 


is 


orientation  of  gages  allowed  for  relatively  easier  compari¬ 
son  between  the  experimental  and  theoretical  values. 

Both  models  were  loaded  incrementally  up  to  a  pressure 
of  500  psi.  Increments  of  50  psi  were  used  for  all  tests. 
Incremental  loading  allowed  time  for  the  scanning  and 
recording  of  the  strain  values  and  pressure  from  the 
measuring  devices.  The  tests  were  all  repeated  at  least 
six  times  in  order  to  minimize  statistical  fluctuations 
and  obtain  consistent  results. 

ROCK  TYPE.  LOCALE  AND  MATERIAL  PROPERTIES 

Two  blocks  of  Precambrian  Milbank  Granite  were 
obtained  from  a  quarry  near  Redfield,  South  Dakota.  A 
detailed  geolgic  description  of  Milbank  granite  is  given  by 
Daniells( 1971) and  White(l973)  .  The  general  geologic 
relationships  and  petrographic  details  of  the  Milbank 
granite  have  been  described  by  Goldich  (19bl)Hall  (1899) and 
Theil  and  Dutton  (1935) 

The  Milbank  granite  is  a  medium-  to  coarse-grained 
light  reddish  brown  rock.  It  consists  of  coarse-grained, 
pink  microcline,  coarse-  to  medium-grained  white  orthoclase, 
medium-grained,  bluish-gray  quartz,  dark  brown  fine¬ 
grained  biotite  and  minor  amounts  of  plagioclase.  A  few 
feldspar  grains  are  found  to  be  larger  than  the  prevalent 


grain  size. 


In  thin  section,  the  rock  is  holocrystalline ,  phaner- 
itic,  medium-  to  coarse-grained  and  inequiangular .  Large 
grains  of  orthoclase  and  microcline  are  characteristically 


surrounded  by  fine-grained,  granulated  quartz  and  ortho¬ 
clase.  The  rock  consists  of  3'/%  anhedral,  irregular,  fine- 
to  medium-grained  quartz,  51%  anhedral,  medium-  to  coarse¬ 
grained  microcline,  10%  anhedral,  coarse-grained  plagio- 
clase,  fine-grained,  irregular  biotite  grains  scattered 
throughout  the  rock  mass,  and  traces  of  magnetite  and 
apatite.  The  minerals  do  not  seem  to  be  preferentially 
oriented;  they  coi  tain  numerous  small  fractures  both  at  and 
across  grain  boundaries. 

Daniells( 1971) has  determined  the  material  properties 
of  Milbank  granite  from  unconfined  compression  tests  on 
1  7/8"  diameter  cylindrical  specimens  in  the  laboratory. 

The  rock  is  isotropic  and  has  an  unconfined  compressive 
strength  of  23,000  psi.  The  modulus  of  elasticity  and 
Poisson's  ratio  of  the  material  are  determined  to  be  8.57  x 
10^  psi  and  0.191  respectively. 

DESCRIPTION  OF  THE  MODELS  AND  EXPERIMENTAL  SET  UP 

The  details  of  the  models  and  the  experimental  set  up 
are  sho^n  in  Figures  18  and  19  Figure  18a  shows  the 


first  model  of  Milbank  granite  of  dimensions  24"x24"x36" 
with  a  cylindrical  hole  of  6  5/8"  diameter  in  the  middle  of 


Schematic  diagram  of  loading  apparatus  and  rock  block  model 


square  face.  The  second  model  has  identical  overall  dimen¬ 
sions  but  the  cylindrical  hole  is  drilled  in  the  middle  of 
one  of  the  rectangular  faces.  These  two  models  were  loaded 
in  uniaxial  compression  by  a  loading  frame  as  shown  in 
Figure  19  This  loading  frame  is  similar  to  the  one  used 
by  Hoskins(  1967 )  for  uniaxial  compressive  loading  of  rock 
and  concrete  blocks. 

The  loading  frame  consisted  of  two  steel  end  plates 
held  at  a  fixed  distance  by  threaded  rods  of  3"  diameter. 
Eight  rods  and  bolts  were  used  to  tie  the  two  end  plates 
together.  The  end  plates  were  1-inch  thick  mild  steel 
plates  with  three  wide  flange  W8x40  beams  24"  long  welded 
to  the  outside  faces  of  the  end  plates.  The  I-beams  were 
stiffened  by  welding  1/2  in.  thick  plates  cut  to  shape  into 
the  webs  at  each  end  of  each  beam.  The  middle  I-beam  at 
one  end  of  the  frame  had  a  7"  pipe  in  the  middle  and  was 
concentric  with  the  6  5/8"  cylindrical  hole  in  the  first 
rock  model.  Large  flat  hydraulic  jacks  (24"x24")  with  and 
without  circular  cut-outs  were  used  to  load  the  rock  and 
these  jacks  acted  against  the  relatively  rigid  frame.  The 
flat  jacks  we  re  sandwiched  between  two  1/4  in.  sh  e  e  t  s  of 
masonite  packing  material  to  obtain  a  uniform  stress 
application  to  the  model. 


MODEL  PREPARATION  AND  TESTING  PROCEDURE 


The  models  were  large  blocks  of  dimensions  24  in.  x 
24  in.  x  36  in.,  and  the  instrumentation  was  installed 
after  fitting  the  model  into  the  loading  frame  as  shown  in 
Figure  19  A  1/4  in.  hole  was  drilled  to  a  depth  of  1/2" 
on  the  borehole  walls  of  the  two  models  at  a  distance  of 
4"  from  the  block  boundary.  In  model  2,  the  hole  was 
drilled  in  a  direction  normal  to  the  loading  axis.  Both 
strain  gage  rosettes  (120°  delta)  and  single  element  gages 
were  used  to  measure  strains  on  the  borehole  walls  and  a- 
round  the  1/4  in.  holes.  The  gages  were  all  1/8  in.  long. 

All  120  degree  rosettes  used  were  Micro  Measurements 
(M-M)  precision  strain  gages,  type  EA-06-125YA-120  with  a 
gage  factor  of  2.065  +1%  and  a  resistance  of  120  ohms  ±0.2%. 
The  single  element  gages  were  M-M  precision  strain  gages, 
type  EA-06-125BT- 120  with  a  gage  factor  of  2.11  +0.5%  and 
a  resistance  of  120  ohms  ±15%.  The  compensating  gage  was 
also  of  this  type.  The  gages  were  a]l  open-faced  general 
purpose  gages  with  <i  polyimide  backing. 

Strain  gage  rosettes  were  installed  on  the  borehole 
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walls  at  locations  diametrically  opposite  to  the  1/4  in. 
holes.  Throe  single-element  gages  were  installed  around 
each  1/4  in.  hole.  These  three  single  gages  constitute  a 
single  strain  rosette.  The  rosette  configuration  chosen 


for  tests  on  models  1  and  2  are  shown  in  Figures  20  and 

It  may  be  seen  that  the  strain  gages  were  all  placed 
at  a  distance  of  1/16"  from  the  hole  boundary.  Since  the 
strain  gradient  around  a  hole  of  this  size  was  compara¬ 
tively  steep,  it  was  necessary  to  devise  a  means  of 
accurately  positioning  the  strain  gages.  Flexible  paper 
templates  were  designed  for  this  purpose.  Both  the  hole 
and  the  positioning  of  strain  gages  were  accurately  drawn 
on  the  template  which  was  then  cut  along  the  hole  circum¬ 
ference  and  gage  boundaries  (Figure  22  ).  The  strain  gages 
were  positioned  in  the  slots  and  then  carefully  attached  to 
the  template  by  taping  along  their  boundaries.  The  lead 
wires  were  next  soldered  on  to  the  gage  terminals  and  the 
entire  unit  was  mounted  on  the  rock  surface. 

The  strain  gages  were  installed  on  the  rock  surface 
using  M-Bond  200,  a  methy 1-2-cyanoacrylate  adhesive  and  M- 
Rond  200  catalyst.  Because  of  its  easy  handling  and 
immediate  room  temperature  cure,  this  adhesive  proved  to  be 
very  satisfactory.  A  coating  of  Dow  Corning  3140  RTV,  a 
clear,  non-corrosive,  flowable,  room  temporciture  curing 
silicone  rubber,  was  applied  over  the  strain  gages  and 
terminal  strips  for  moisture  resistance  and  to  avoid  gage 
instability  and  any  mechanical  damage. 

Before  bonding  the  strain  gage  on  the  borehole  wall 


Fig.  21  Trepanning  hole  rose the  configuration  in  model  ? 


it  was  necessary  to  prepare  the  rock  surface.  The  surface 
was  first  sanded  and  then  thoroughly  cleaned  with  acetone. 

A  generous  spray  of  Freon  degreaser  followed.  The  rosette 
gages  were  then  bonded  to  the  surface  with  the  M-Bond  200 
adhesive  and  catalyst.  Excessive  adhesive  was  forced  from 
under  the  gages  by  pressing  on  the  gages  from  the  center 
outward.  This  reduced  the  possibility  of  air  bubbles  being 
present  under  the  gages.  Immediately  upon  completion  of 
wipe-out  of  the  adhesive,  firm  thumb  pressure  was  applied 
to  the  gages  and  terminal  area  for  several  minutes  after 
which  the  gages  were  solidly  bonded  in  p] ace .  The  gages 
were  then  visually  chocked  for  improper  bonding.  Gage-to- 
gage  specimen  resistance  war  measured  with  a  multimeter  to 
make  sure  the  gages  were  intact. 

An  identical  gage  bonded  in  the  same  way  to  a  separate 
block  of  material  was  used  as  the  temperature  compensating 
gage.  Only  one  compensating  gage  was  used  throughout  the 
experimental  work. 

Two  types  of  strain  .indicators  were  used  to  record 
strain  readings.  A  Bean  Model  201  Digital  Ft  rain  Indicator 
in  association  with  a  5-channel  Model  301  switch  and  balanc¬ 
ing  unit  were  used.  Both  indicators  proved  stable  and 
reproducibility  of  strain  readings  was  good. 

Two  flat  jacks  were  used  to  apply  normal 


stress  to  the  models. 


Pressure  was  applied  to  each  jack  by 
a  3000  psi  hand  pump.  The  pressure  applied  by  the  hand 
pump  was  monitored  by  the  previously  calibrated  5000  psi 
pressure  gage. 

The  lead  wires  from  the  strain  gages  were  connected 
to  the  strain  indicators.  The  indicators  were  then  connected 
to  power  supplies  and  allowed  to  warm  up  for  a  period  of 
severa.  ..linutes.  During  the  warm  up  time,  all  the  connec¬ 
tions  were  checked.  Another  check  on  the  pressure  gage  and 
pressure  valve  was  made. 

The  initial  strain  values  of  all  strain  gages  were 
checked  and  recorded.  Pressure  was  then  applied  through 
the  hand  pump  in  50  psi  increments.  At  each  increment,  the 
strain  values  of  all  strain  gages  were  recorded.  The  process 
was  continued  until  a  pressure  of  500  psi  was  applied. 

Then  the  pressure  was  released  gradually  until  it  was 
brought  back  to  zero.  The  strain  values  were  again  recorded 
at  every  50  psi  pressure  drop  during  unloading  process.  The 
loading  and  unloading  cycle  was  repeated  at  least  six  times 
and  the  final  strain  values  are  calculated  by  averaging 
the  12  sets  of  readings  so  obtained. 

COMPARISON  AND  D IS CUSS ION  OF  RUSULTS 

The  average  strain  values  obtained  from  models  1  and 
2  are  presented  in  Tables  4  and  5  The  difference 
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between  the  strains  around  the  hole  and  tno  corresponding 
strains  measured  at  a  diametrically  opposite  point  gives 
the  relaxation  strains  at  these  points.  Therefore  the 
1/4  in.  hole  is  termed  a  trepanning  hole  although  no  rosette 
was  undercored.  This  procedure  also  eliminates  any  residual 
stresses  present  in  the  rock  and  allows  a  direct  comparison 
to  be  made  between  the  experimental  results  and  predicted 
va lues . 

From  Table  5  it  may  be  seen  that  the  response  of 
the  gages  is  linear.  The  values  in  the  first  three  columns 
are  tne  strains  resulting  from  a  homogeneous  compressive 
stress  field.  For  an  applied  stress  (o  )  of  500  psi  along 
the  borehole  axis,  the  principal  stresses  and  strains 
computed  from  experimental  values  using  elastic  stress- 
strain  relationships  are  compared  (Table  7  )  with 
Hiramatsu's  elasticity  solution  (1962)and  three-dimensional 
finite-element  solution.  Doth  theoretical  so]  itions  give 
identical  results  which  is,  of  course,  expected  for  this 
simple  uniaxial  compressive  field.  The  experimental  values 
show  good  agreement  in  general.  However,  the  discrepancy 
is  as  large  as  15%  for  minimum  stress  and  strain  values. 

This  can  be  partly  accounted  for  by  the  fact  that  relatively 
small  strains  were  measured  and  the  strain  .indicator  readings 
can  only  be  estimated  within  2.5  nicroinches/j  nch .  Therefore 
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any  small  discrepancy  in  the  measured  strains  will  result 
in  a  very  large  percentage  for  smaller  measured  strains. 

TABLE  7 

COMPARISON  BETWEEN  EXPERIMENTAL  AND  THEORETICAL  SOLUTIONS 
(MODEL  1,  MEASUREMENTS  ON  SIDEWALL),  p=500  psi 


Principal  Stresses 
and  Strains 

Experimental 

Hiramatsu ' s 
Solution 

3D-FEM 

a 

+  4 . 0 

0 

0 

max 

o  . 

-427.74 

-500.0 

-500.0 

min 

e 

max 

-10 .00 

-11.14 

-11.14 

e  . 
mm 

-50.0 

-58.34 

-58.34 

The  strains  measured  around  the  trepanning  hole  in 
model  1  are  given  in  columns  4  to  6  of  Table  5  By 
subtracting  the  corresponding  strains  in  columns  1  to  3 
from  these  values,  the  relaxation  strains  were  obtained  and 
are  shown  in  columns  7  to  9 .  The  relaxation  strains  at  the 
three  gage  points  were  also  predicted  by  performing  a  three- 
dimensional  substructure  finite-element  analysis  on  model  1. 
The  theoretical  strains  were  computed  by  taking  weighted 
averages  over  the  elements  on  which  the  gages  were  placed. 
Predicted  strain  values  versus  measured  strain  values  are 
shown  in  Figure  23  it  can  be  seen  that  the  experimental 
values  lie  very  close  to  the  predicted  values  for  gage  1. 

The  agreement  is  quite  good  for  gages  2  and  3  but  the 
percentage  discrepancy  is  large  at  some  discrete  points. 
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Fig.  23  Comparison  between  experimental  and 


predicted  relaxation  strains  for  model  1 
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experimental  values. 

From  the  experimental  relaxation  strain  values,  the 
applied  stress  field  was  calculated  using  Russell's  two- 
dimensional  data  reduction  computer  program  for  deep  hole 
device.,  The  stresses  along  and  normal  to  the  bore¬ 

hole  axis  were  found  to  be  —487  psi  and  -27  psi  respective¬ 
ly.  To  obtain  the  applied  stress  field  from  the  three- 
dimensional  analysis,  a  total  of  nine  relaxation  strain 
measurements  were  required  whereas  only  throe  measurements 
were  xaade .  Since  the  model  is  subjected  to  a  uniform 
homogeneous  compressive  stress  field,  the  same  three 
relaxation  strains  should  exist  at  each  and  every  G-position, 
if  there  are  no  statistical  fluctuations  in  the  experimental 
values.  Using  equation  (3.24b)  with  a  unit  weight  matrix 
and  equation  (3.31),  the  principal  stresses  were  found  to 
be  equal  to  -486  psi,  -23  psi  and  -23  psi .  The  agreement 
is  excellent  since  the  difference  in  the  applied  stress 
field  (oz=500  psi)  is  only  2.4?4. 

In  testing  model  2,  the  linear  response  of  the  model 
in  the  designed  loading  range  was  first  verified  by  the 
strains  measured  on  the  sidewall  of  the  borehole.  Figure 
24  shows  the  response  of  strain  rosette  to  the  uniaxial 
compressive  stress  applied  normal  to  the  borehole  axis. 

The  dispersion  of  strains  for  any  single?  gage  clearly 
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As  already  mentioned,  this  is  mainly  because  of  the  very 
small  strain  values. 

In  addition  to  discrepancies  due  to  differences  in 
recording  smaller  strains,  discrepancies  could  have  been 
introduced  duo  to  any  of  the  following  causes: 

(1)  Inaccuracies  in  determining  the  elastic  constants. 
The  predicted  strains  are  dependent  on  the  Young's  modulus 
and  Poisson's  ratio  of  the  material.  So,  any  change  in 
these  values  directly  affects  the  strains. 

(2)  The  measurements  were  made  at  points  not  very  far 
from  the  outside  boundary  of  the  block  of  rock.  This  was 
necessary  to  make  the  drilling  process  simple  and  accurate, 
but  the  ideal  location  would  bo  within  the  middle  third 
portion  of  the  model. 

(3)  The  rock  satisfies  the  homogeneity  condition  only 
in  a  statistical  sense.  Any  of  the  gages  may  cover  a 

particular  crystal  whose  elastic  constants  differ  from  those 
of  the  model  material. 


(4)  Another  possible  source  of  error  may  be  due  to  the 
inadequacy  of  the  finite-element;  model.  Because  of  the 
limitation  imposed  by  the  computer  core  storage,  models 
with  a  much  finer  subdivision  were  not  investigated.  A 


finer  mesh  or  a  model  with  refined 


(hi  iher;  order)  elements 


might  result  in  a  much  closer 


agreement  with  the 
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indicates  a  linear  trend. 

From  the  strain  values  in  the  first  three  columns  of 
Table  6  the  principal  stresses,  strains  and  their 
directions  were  determined  using  the  equations  derived  in 
Appendix  B.  If  gage  1  was  made  to  coincide  with  the  direc¬ 
tion  of  the  borehole  axis,  the  maximum  principal  strain 
direction  would  be  0°  with  reference  to  gage  1.  But  the 
gage  was  inclined  at  approximately  15°  to  the  borehole  axis, 
and  with  an  applied  compressive  stress  o  =-500  psi,  the 
maximum  principal  strain  direction  was  found  to  be  equal  to 
-12°53 1  (clockwise  negative).  The  principal  strain  direc¬ 
tion  was  determined  for  every  set  of  strain  readings  and 
the  results  were  plotted  as  shown  in  Figure  25  From  this 
figure,  it  may  be  seen  that  the  principal  strain  direction 
does  not  show  a  variation  of  more  than  +2°  which  proves 
that  the  strains  were  highly  consistent  and  responded 
linearly  to  the  applied  stress  field. 

For  an  applied  compressive  stress  field  of  500  psi, 
the  principal  stresses  and  strains  on  the  sidewalls  of  a 
borehole  were  computed  using  Howland's  solution  (1930) 
Hiramatsu's  three-dimensional  elasticity  solution  (1962)nnc] 
three-dimensional  finite-element  solution.  In  Table  8 
these  results  are  presented  along  with  experimentally 
determined  values.  Howland's  two-dimensional  elasticity 
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solution  for  a  semi-infinite  strip. 

TABLE  8 


COMPARISON  BETWEEN  EXPERIMENTAL  AND  THEORETICAL  SOLUTIONS 


(Model 

2,  Measurements 

on  Sidewall,  p=-500 

psi) 

Principa 1 
Stresses 

and  Strains 

Experimental 

Howland 

Hiramatsu 

3D-FEM 

a  .  (psi) 
min 

-1767.49 

-1660.0 

-1500. 

-1664.0 

o  (psi) 

max 

-20.45 

- 

0 

-5.21 

a  .  (psi) 
nun 

-205.79 

-193.70 

-179.68 

-194.72 

o  (psi) 

max 

+37.01 

- 

+21.70 

+  36.5 

under  tension  and  the  three-dimensional  finite- element 
solutions  show  the  closest  agreement  with  the  experimentally 
determined  values.  The  maximum  principal  stress  and  strain 
values  cannot  be  obtained  from  Howland’s  solution  which  is 
two  dimensional.  Iliramatsu's  solution  which  is  valid  only 
if  the  block  may  be  reyarded  as  indefinitely  extended  in 
two  dimensions,  yjves  a  principal  stress  value  of  -1500  psi 
corresponding  to  a  stress  concentration  factor  of  3.0.  Both 
Howland’s  solution  and  finite-element  solution  yive  a  some¬ 
what  higher  stress  concentration  factor  of  3.3  at  this 
point.  The  three-dimensional  finite-element  solution 
shows  excellent  agreement  with  the  experimentally  determined 
values,  the  difference  being  less  than  6%  in  all  but  one 


case . 
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The  strains  around  the  trepanning  hole  are  shown  in 
columns  4  to  6  of  Table  6  For  obtaining  the  relaxation 
strains  for  model  2,  the  strains  along  vertical  and  +120° 
directions  had  to  bo  computed  because  of  the  positioning 
of  the  rosette  at  this  location.  The  relaxation  strains 
thus  obtained  from  experimental  values  were  plotted  against 
the  values  predicted  by  the  three-dimensional  substructure 
finite-element  analysis  in  Figure  26  The  theoretical 
values  wore  once  again  determined  by  taking  weighted 
averages  of  the  strains  in  elements  over  which  the  gages  were 
placed.  Very  good  agreement  has  been  found  and  the  differ¬ 
ences  seldom  exceed  +_15%. 

CONCLUDING  REMARKS 

The  results  of  the  experimental  investigation  have 
confirmed  the  application  of  three-dimensional  substructure 
analysis  presented  in  the  previous  chapter  for  determining 
primary  stresses  in  rock  from  measurements  on  the  sidewalls 
of  boreholes.  This  has  been  demonstrated  by  the  excellent 
agreement  between  the  experimental  and  theoretical  values 
of  the  relaxation  strains  in  the  tested  models.  With  the 
experimental  set  up  that  was  used  and  the  magnitudes  of 
strains  that  were  measured,  agreement  between  theoretical 
and  experimental  values  within  +10  to  15%  would  be  consider¬ 
ed  acceptable.  From  the  values  presented  in  Tables  ? 
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and  8  and  Fiyures  23  and  26,  it  may  be  pointed  out  that 
the  differences  between  experimental  and  theoretical  values 
are  less  than  +10%  in  a  majority  of  cases. 
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ANALYTICAL  METHOD 

The  purpose  of  this  section  is  to  describe  the  results 
of  an  investigation  into  a  suitable  means  of  converting  the 
nine  components  of  relieved  strain  into  principal  stresses. 

These  stresses  are  assumed  to  act  a  point  in  the  center  of 
the  hole  at  the  depth  of  the  device  prior  to  drilling  the  bore¬ 
hole  and  the  trepanning  holes.  A  further  objective  of  the 
analysis  is  to  provide  information  useful  in  designing  an 
efficient  device.  It  is  assumed  that  the  mechanical  properties 
of  the  rock  in  the  region  of  the  measurement  will  be  known  from 
either  laboratory  tests  on  core  cut  from  the  borehole  or  from 
a  device  capable  of  measuring  mechanical  properties  in  place. 

Assumptions  In  order  to  estimate  the  stress  field  at 
a  point  that  existed  prior  to  drilling  the  borehole  and  the 
trepanning  holes,  a  number  of  assumptions  are  made.  These 
assumptions  may  be  divided  into  those  inherent  in  the  strain 
relief  system  and  those  made  to  simplify  the  analysis  of  the  data. 

Assumptions  inherent  in  the  strain  relief  method  and 
device  are  1)  strain  changes  resulting  from  residual  stresses 
are  negligible  in  comparison  with  those  due  to  gravity,  active 
tectonic  forces  and  any  other  forces  presently  acting,  (Residual 
stresses  are  defined  as  those  acting  in  a  body  ol  rock  in  the 
absence  of  any  external  force  or  force  field),  2)  stress  gradients 
are  small  enough  so  that  they  cause  only  negligible  chances  in 
strain  over  the  dimensions  used  in  the  device,  in  other  words,  the 
difference  in  stress  between  the  highest  and  lowest  strain  gage 
units  is  negligible  in  comparison  with  the  total  stress,  3)  strain 
gages  bond  to  continuous  rock  and  do  not  span  significant  voids 
or  joints  and  4)  the  device  anchoring  jacks,  bonding  of  the  strain 
gages  and  drilling  of  the  holes  do  not  significantly  alter  the 
behavior  of  the  rock  in  the  neighborhood  of  the  measurement,  i.e., 
tests  for  mechanical  properties  run  on  cores  should  be  representative 
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of  the  in  situ  rock. 

Assumptions  made  to  simplify  the  analysis  arc  1)  the 
rock  in  the  neighborhood  of  the  trepanning  holes  may  be 
modeled  as  a  linear  isotropic,  homogeneous  elastic  materi¬ 
al,  2)  any  inelastic  effects  that  may  occur  very  close  to 
the  holes  are  not  significant,  3)  the  trepanning  holes  act 
independently,  i.o.,  drilling  the  first  trepanning  hole 
does  not  affect  the  results  from  drilling  the  second  and 
third  holes,  4)  the  curvature  of  the  borehole  may  be  neg¬ 
lected  in  computing  the  principal  stresses  from  the 
relieved  strains  and  5)  the  direction  parallel  to  the  bore¬ 
hole  is  a  principal  direction.  Assumption  number  4, 
neglect  of  borehole  curvature  in  comparison  with  tre¬ 
panning  hole  diameter,  allows  us  to  use  the  two-dimensional 
plane  stress  solution  lor  determing  stresses  from  relieved 
strains . 

Method  Using  the  above  assumptions,  an  approximate 
solution  may  be  found  by  the  following  procedure: 

1)  Analyze  the  striin  relief  for  each  strain  gage  unit  by 
using  a  modification  of  a  method  published  by  Soete  (1950). 

2)  Check  the  results  from  the  three  trepanning  holes  for 
consistency  in  the  vertical  stress  in  accordance  with 
assumption  5  above,  if  assumption  5  is  satisfied,  proceed 
with  the  analysis,  otherwise,  the  approximate  analysis 
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O 

f  a i Is, 

3)  Since  the  direction  of  the  borehole  is  a  principal 

; 

direction,  the  second  and  third  principal  directions  at 
each  trepanning  hole  must  lie  in  a  plane  normal  to  the 
axis  of  the  borehole,  therefore,  we  have  the  tangential 
component  of  si i ess  at  three  points  around  the  circumfer¬ 
ence  of  the  borehole  spiced  at  120  degrees  and  assumed  to 
bo  acting  in  Liie  same  plane, 

4)  Since  the  vertical  stress  (assumed  not  to  vary  with  X 
and  Y)  can  lie  shown  to  have  no  effect  on  the  horizontal 
stress  distribution,  the  three  tangential  stresses  at  120 
degrees  are  sufficient  to  determine  the  two  principal 
stresses  and  the  principal  direction  in  the  horizontal 
plane.  In  tl  i  next  two  sections,  the  modified  Soete  method 
and  the  principal  stresses  in  the  horizontal  plane  are 
developed . 

Strain  lb- Inf  Near  a  Trepanning  Hole  The  stress  and 
strain  fields  e: isting  in  a  material  will  be  disturbed  by 
the  drilling  of  a  trepanning  I  ole.  The  change  in  the 
strain  field  n<  ,ir  the  hole  can  be  measured  by  means  of 
electrical  resistance  strain  gages  or  some  other  suitable 
means.  The  basic  relationship  between  the  original  strain 
field,  the  change  in  strain  and  the  final  strain  field  may 
be  expressed  by  the  following  equation. 
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Final  Strain  Field  =  Original  Field  +  Change  in 
Field  (1) 

Since  strain  gages  measure  only  an  average  strain  over 
their  length  and  not  the  entire  strain  field,  we  can 
rearrange  Eq  .  1  and  interpret  it  to  hold  at  any  point  in 
the  strain  field.  Then  for  any  point, 

Measured  Strain  =  Final  Strain  -  Original  Strain  (2) 
where  measured  strain  has  been  equated  with  the  average 
change  in  strain  over  the  length  of  the  strain  gage  used. 

If  we  assume  that  stress  gradients  will  produce  only 
negligible  changes  in  stress  over  the  dimensions  involved, 
the  2-D  stat<  of  stress  in  the  borehole  wall  (assumed  to 
be  plane)  may  be  expressed  in  terms  of  two  principal 
stresses  and  a  principal  angle.  Similarly,  the  final  state 
of  stress  around  the  trepanning  hole  may  be  expressed  in 
terms  of  the  same  principal  stresses  and  angle  by  using 
the  2-D  plane-stress  solution  originally  developed  by 
Kirsch  (1898).  Since  we  arc  neglecting  the  curvature  of 
the  borehole  for  the  analysis  of  the  strain  relief  at  the 
trepanning  holes,  we  are  essentially  on  the  surface  of  a 
half-space  where  the  normal  component  of  stress  is  zero. 
Consequently,  the  plane  stress  solution  is  the  appropriate 
approximation . 


The  above  argument  implies  that  the  measured  strain 
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at  a  point  near  a  trepanning  hole  may  be  expressed  as  a 
function  of  two  principal  stresses  and  a  principal  angle 
acting  in  a  plane  normal  to  the  axis  of  the  trepanning  hole. 
Thus 

Measured  Strain  =  Function  (J|,  u)  (3) 

where  o,  and  are  the  principal  stresses  and  «  is  the 
principal  angle  all  acting  in  a  plane  normal  to  the  axis  of 
the  trepanning  hole  being  considered.  As  shown  below  the 
function  implied  by  Fq.  3  is  nonlinear.  Therefore,  if  we 
apply  Fq.  3  to  each  of  the  three  strain  gages  around  a 
trepanning  hole,  we  will  have  a  set  of  throe  simultaneous, 
nonlinear,  al  jebraic  equations  which  may  be  used  to  solve 
for  the  three  unknowns  5 , ,  o ;  and  0  acting  near  the  tre¬ 
panning  hole. 

In  equation  form,  from  Socte  (1950),  we  have 


L  *  =  A 

o  J 

V 

i  +  h 

V.  Vj 

(4a) 

=  A 

f  i 

v,  + 

L  V.cos  2(9+0 

i  t 

(4b) 

c  '  =  A 

-  U 

V  i  * 

13  V,  cos  2  ( o-  i) 

-  t 

(4c) 

where  <■ 

i 

s  the 

angle  between  gages  as  shown  on  Fig.  3, 

V, 

= 

(a  !  + 

.  )/r: 

(5a) 

V 

= 

(j  i" 

, .  )  /b 

(5b) 

V3  =  cos(2H)  (5c) 

L  =  Young's  modulus  of  elasticity 

e "  =  avei  tgo  measured  strain 
m 

1  •  r 

*  — — —  ng  t  *  dr  (used  to  develop  f.qs.  (4)  &  (G), 

in?  i;ii  L 

nu 
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where  m  =  0,  +u,  and  -a, 

rm.  an<3  are  respectively  the  radial  distances  from  the 

center  of  the  hole  to  the  inner  and  outer  edges  of  the 
strain  gage  grid  at  angle  m, 

c"  is  the  function  representing  the  variation  of  the 
radial  strain  with  the  radial  coordinate,  r,  at  the  angle 
m, 


A 


m 


1+y  a^ 

Err 

nil  nu 


B 

m 


mi  ni _ 


1-1 


+  (1+v) 
4 


a ‘ ( r ‘  +r 


r  +r '  ) 
m  i  m  i  m ?  in2 


r /  r  • 
m  t  m  2 


(6a) 

(6b) 


m  =  0,  +u,  -a,  and 

a  =  radius  of  trepanning  hole. 

Note  that  the  function  representing  the  variation  of  radial 
strain  is  averaged  only  over  the  radial  direction,  this 
implies  that  the  width  of  the  strain  gage  grid  must  be 
small  in  comparison  with  the  radius  of  the  trepanning  hole. 

The  three  unknowns  in  Eqs.  4  are  V  ,  V  ,  and  V  .  If 


1  2  3 

the  three  values  of  V  are  known,  we  can  easily  solve  for 
0  -and  6.  Before  proceding  with  the  solution,  the 

cosine  terms  in  Eqs.  4b  and  4c  must  be  expressed  in  terms 
of  V]  =  cos  2«.  This  can  be  accomplished  by  using  standard 
trigonometric  identities.  Then 


AV  +  B  V  V  =  e  ' 
O  i  0  2  3  0 


(7a) 


A  V  +B  V  (V  cos  2U-/1-V2  sin  2a)=e; 

+  Ct  i  —a  2  3  3  +a 


(7b) 


a 
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A  V  +B  V  (V  cos  2a+v  1-V  7  sin  2a)=e'  (7c) 

-  a  i  -  c*  ?  3  3  -  x 

Since  the  right  hand  sides  of  Eqs.  7  are  the  neasured 
changes  in  strain  at  the  three  angles  and  all  of  the  sym¬ 
bols  on  the  left  hand  side,  except  the  three  V's,  can  be 
computed  from  a  knowledge  of  the  material  properties  and 
the  geometry,  Eqs.  7  constitute  a  set  of  three  simultaneous, 
nonlinear,  algebraic  equations  in  terms  of  the  three 
unknown  values  of  V. 

Eqs.  7  differ  from  those  presented  by  Soete  (1950). 
Soete  assumed  that  each  strain  gage  was  at  exactly  the  same 
distance  from  the  center  of  the  trepanning  hole.  In  that 
case,  the  three  A's  in  Eqs.  7  could  be  represented  by  a 
single  symbol  A'.  Similarly,  the  three  B's  in  Eqs.  7 
could  be  represented  by  B'.  The  simpler  geometry  allows 
the  direct  solution  of  the  equations  and  is  given  by  ^ete. 
As  v/ill  be  shown  later,  the  principal  stresses  are  very 
sensitive  to  small  changes  in  geometry.  Therefore,  a 
solution  was  constructed  for  the  more  general  Eqs.  7. 

To  construct  a  solution  to  Eqs.  7,  we  first  add  Eq . 

7b  and  Lq .  7c  to  eliminate  the  square  root  term.  The  pro¬ 
duct  V 2^ 2  remains  but  can  be  eliminated  by  using  Eq.  7a  to 
give 

c  +t  -Be 

V  — ±£i — — —  (8a) 

1  A  +A  -BA 
-  v  -a  o 

where  B  =  (B1  +B  )cos(2a)/B 
+  a  -or  o 

An  expression  may  be  obtained  for  by  first  sub¬ 
tracting  Eq.  7c  from  Eq.  7b  and  then  eliminating 
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V  by  using  Eg.  7a.  After  some  algebraic  manipulation, 

3 

(A  -A  )V  -it'  -e'  )2  2  . 

V  =  { — — - — - 1 - - - -  +(l'-A  V  )  }  */B  (8b) 

2  B  tan  2a.  °  °  1  ° 

The  remaining  unknown  may  be  found  from  Eg.  7a  to  be 

V  =  (c  '-A  V  )  /  ( B  V  )  (8c) 

3  O  l  0  2 

It  should  be  noted  that  Eqs.  8  must  be  solved  sequentially 
since  the  expression  for  V2  involves  and  the  expression 
for  involves  both  and  V  ,y .  Algebraic  substitutions 
could  be  made  but  the  sequential  calculation  causes  no 
difficulty  and  appears  to  be  the  most  efficient  way  of 
outaining  numerical  values  for  the  unknowns. 

After  calculating  the  values  for  V^,  \’0  and  V^,  Eqs.  5 
are  used  to  find 


c  =  E(V.+V,)/2 
i  I  ~ 

(9a) 

0  =  E(V1-V2)/2 

(9b) 

,  -1 

(9c) 

o  =  +1/2  cos  (V^) 

The  correct  sign  for  the  principal  angle  9  may  be  found  by 
substituting  into  either  Eq .  4b  or  4c.  If  the  equation  is 
satisfied,  the  correct  sign  has  been  chosen  for  0.  The 
angle  0  is  defined  as  the  clockwise  angle  measured  from  the 
direction  of  principal  stress  o,  to  the  direction  chosen  as 
0.  It  should  be  noted  that  the  equations  derived  above 
assume  that  tension  is  positive.  Also,  the  a  derived  in 
the  above  equations  may  not  be  the  algebraic  maximum 
principal  stress. 

Principal  Stresses  in  3b.  The  strain  change  data  from  each 
of  the  three  trepanning  holes  are  analyzed  using  Eqs.  8  and 


and  9.  As  previously  stated,  if  the  results  from  each 
trepanning  hole  do  not  indicate  that  the  vertical  direction 
is  principal,  the  approximate  method  fails  and  a  more  exact 
approach  must  be  taken,  see  Sateesha  (1974).  It  is  anti¬ 
cipated  that  in  the  majority  of  practical  situations,  the 
vertical  direction  will  be  vertical  and  the  average  of  the 
three  vertical  stress  components  from  the  three  trepanning 
holes  will  lie  an  adequate  representation  of  the  magnitude 
of  the  vertical  principal  stress.  At  this  stage  we  know 
one  principal  direction  „s  vertical  and  it  follows  that  the 
remaining  two  principal  directions  are  in  the  horizontal 
plane.  It  remains  for  us  to  solve  for  the  two  principal 
stress  magnitudes  and  one  principal  angle  in  the  horizontal 
plane.  From  the  analysis  of  the  three  trepanning  holes, 
tangential  stresses  are  now  known  at  throe  points  around 
the  circumference  of  the  borehole  as  shown  on  Fig.  27 

Again  the  problem  may  be  solved  by  using  the  Kirsch 
solution 

Jt  =  (P1+P2)_2 (p]"p2)  COS  23  (10) 

where  u  is  a  t  ingential  component  of  stress  at  the  bore¬ 
hole  surface,  p^  and  p  are  principal  stresses  acting  in 
the  horizontal  plane  prior  to  drilling  (tension  positive) 
and  3  is  the  principal  angle  measured  clockwise  from  the 
direction  of  p  t  o  tire  reference  direction  0.  I  f  we  denote 
the  three  known  tangential  stress  components  at  the  bore¬ 
hole  surface  as  a  where  y  is  generally  120 

degrees,  and  apply  the  Kirsch  solution  three  times,  we  have 
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II 

o 

o 

'  V2U2U3 

(11a) 

Jt+y 

=  U.. -2L. cos2 (8+y) 

(lib) 

°t-y 

where 

=  U1-2U2cos2(8-y) 

(11c) 

ui  - 

Pl  +  P2 

(12a) 

c 

ro 

II 

P 1  ”  P  2 

(12b) 

U3  = 

cos  2B 

U2c) 

Eqs.  (11)  have  exactly  the  same  form  as  Eqs.  (4)  if  we  make 
the  following  associations 


U.-V. 
l  l 


i  =  1,  2,  3 


a.  •+  £ 
to  o 


l-'A 


m 


Jt+y"fc+a 


-2  *B 


m 


a.  ►£ 

t- y  -a 


y-^a 


Using  Eqs.  8  and  die  above  associations,  the  solution  to 
Eqs.  (11)  is 


U  = 


U  = 

2 

u3  = 


o .  ,  +o .  -2  cos  2y 

t+f  t-y _ _ 

2 (1-cos  2y 


o.  -o^  2 

-< < .  t-i,  t-r  {  + 
2  sin  2 y 

(U  -u  )/(2U  ) 

1  tO  2 


(  0  -u  )  }  V2 
to  1 


(13a) 

(13b) 

(13c) 


rrom  Eqs. 

(12) ,  the  principal 

stresses  and  angle  in  the 

horizontal 

p’ane  are 

P1  “ 

(U1+U2)/2 

(14a) 

P2  = 

(U1-l,2)/2 

(14b) 

8  = 

- h  cos  i(V3) 

(14c) 

J 
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The  appropriate  sign  for  0  may  be  found  by  substituting  0 
with  an  assumed  sign  into  Eq.  (lib)  or  Eq .  (11c),  if  the 

correct  sign  is  assumed,  the  equation  is  satisfied. 

RESULTS 

l°n  of  Results  The  above  equations  for  determin¬ 
ing  principal  stresses  can  be  checked  for  algebraic 
consistency  by  assuming  a  stress  field,  computing  the 
corresponding  changes  in  strain  from  Eqs .  4  and  then 
proceeding  through  the  analysis  as  though  these  were  actual 
measured  strains.  If  the  equations  are  consistent,  the 
original  stress  field  should  be  recovered.  This  algebraic 
check  has  been  performed  but  will  not  be  presented  here  in 
order  to  conserve  space. 

The  algebraic  check  gives  no  information  on  the 
validity  of  the  assumptions  inherent  in  the  method.  These 
assumptions  can  be  checked  by  either  performing  a  more 
exact  analysis  and  comparing  results  or  by  attempting  to 
verify  the  results  experimentally.  Sateesha  (1974)  has 
performed  a  limited  experimental  study  in  the  laboratory. 

A  block  of  Milbank  Granite  2  ft.  x  2  ft.  in  cross- 
section  and  i  ft.  long  was  drilled  leaving  a  6  inch 
diameter  borehole  parallel  to  the  3  ft.  dimension.  The 
block  was  placed  in  a  specially  designed  loading  frame  and 
loaded  parallel  to  the  borehole  with  a  pressure  of  500  psi . 


Strain  gages  placed  on  the  borehole  wall  approximately  6 
inches  from  the  end  registered  the  longitudinal  and  tangen¬ 
tial  components  of  strain.  The  block  was  then  unloaded  and 
a  *4  inch  diameter  trepanning  hole  was  drilled  anto  the 
borehole  wall.  The  block  was  again  loaded  and  the  strains 
around  the  trepanning  hole  were  recorded.  The  difference 
between  the  strain  readings  without  the  trepanning  hole 
and  those  with  the  hole  were  taken  to  be  the  change  in 
strain  that  would  have  been  recorded  if  the  trepanning  hole 
had  been  drilled  while  the  block  was  under  load.  The 
advantage  of  the  above  experimental  procedure  is  that 
residual  stresses  arc  eliminated  from  the  strain  changes. 
This  is  important  since  the  strain  level  is  relatively  low 
due  to  the  low  stress  level  and  the  relatively  high  modulus 
of  elasticity  of  the  granite  used. 

Results  obtained  from  the  test  are  shown  on  Tabic  9 
The  agreement  shown  is  relatively  good  and  the  discrepan¬ 
cies  are  well  within  the  tolerances  of  the  strain 
measurements  and  the  material  properties. 


Applied  Stress 
(psi ) 


Stress  Computed  from  Measured 
Strain  Changes 
p-i  ) 


o  500  487 

0  0  27 

x 

0  0  27 

y 

Table  9  Comparison  of  applied  stresses  with  those 
computed  from  measured  strain  changes. 

The  results  from  one  laboratory  test  certainly  do  not 

verify  all  the  assumptions  inherent  in  the  technique  itself 

or  the  approximate  method  of  analysis.  For  example,  the 

question  of  what  is  the  smallest  ratio  of  borehole  diameter 

to  trepanning  hole  diameter  for  which  the  neglect  of  the 

borehole  curvature  is  appropriate  has  not  been  answered. 

Nevertheless,  the  results  are  sufficiently  encouraging  to 

warrant  further  study  of  the  equations  to  determine  the 

sensitivity  of  the  results  to  errors  in  Poisson's  ratio  and 

gage  pattern  geometry. 

Sensitivity  and  Error  Analysis 

As  noted  in  the  introduction,  one  of  the  objectives  of 

this  study  is  to  provide  meaningful  input  to  the  design 

process.  One  of  the  basic  questions  that  arises  j.n  the 

design  of  the  device  concerns  the  most  efficient  size  of 

trepanning  hole  that  is  practical  and  what  size  of  strain 

gage  should  be  used  in  conjunction  with  the  trepanning 
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hole.  If  the  problem  is  viewed  from  the  instrumentation 
point  of  view,  it  is  apparent  that  a  relatively  large 
signal  is  desired  in  order  to  minimize  the  inevitable  noise 
problems  in  transmitting  the  signal  from  the  device  to  a 
recorder  at  the  surface.  Consequently,  an  efficiency 
appropriate  to  this  situation  can  be  defined  as  the  ratio 
of  strain  relieved  in  a  particular  direction  (usually 
vertical)  by  trepanning  to  the  maximum  strain  that  could 
occur  in  that  direction  in  a  uniaxial  test,  a/E.  Taking 
a  vertical  stress  of  1000  psi  and  a  horizontal  stress  of 
250  psi  and  maintaining  a  constant  distance  between  the 
edge  of  the  hole  and  the  inner  edge  of  the  strain  gage  grid 
°f  1/16  inch,  relieved  strains  have  been  calculated  for 
various  hole  diameters  and  strain  gage  grid  lengths  using 
Eqs.  4.  The  results  of  Efficiency  vs.  Grid  Size/Hole 
Diameter  are  shown  on  Fig.  28  Fig.  28  indicates  that  for 
maximum  efficiency  the  ratio  of  grid  size  to  hole  diameter 
should  be  as  small  as  possible.  This  is  reasonable  since 
a  strain  gage  averages  the  strain  over  its  length  and  the 
distribution  of  radial  strain  is  nonlinear  with  the  highest 
values  occurring  adjacent  to  the  hole.  The  efficiencies 
shown  hold  only  for  the  value  of  Poisson's  ratio  used  and 
the  fixed  clearance  between  the  strain  gage  grid  and  the 
edge  of  the  hole.  Nevertheless,  the  conclusion  seems  to 


Efficiency 
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be  warranted  and  the  values  approximate  for  the  practical 
situation.  A  ratio  of  grid  size  to  hole  diameter  of  1/6 
would  coincide  to  a  1/16  inch  grid  and  a  hole  diameter  of 
6/16  =  3/8  inch.  For  a  nominal  borehole  of  6  inch  diameter, 
these  values  of  trepanninn  hole  diameter  and  strain  gage 
grid  size  are  perfectly  reasonable  although  practical 
constraints  on  the  friction  bonding  strain  gages  may  pre¬ 
clude  the  use  of  these  relatively  small  gages. 

Another  question  worthy  of  consideration  relates  to 
how  sensitive  the  computed  stresses  are  to  errors  in  the 
strain  gage  pattern  radius.  Assuming  the  same  stress  field 
as  noted  above,  the  per  cent  error  in  vertical  and  hori¬ 
zontal  components  of  stress  may  be  computed  for  the  assumed 
correct  position  of  the  strain  gages  and  for  gage  pattern 
radii  differing  from  the  assumed  by  various  amounts.  The 
results  are  presented  on  Fig.  29  which  shows  the  relative 
error  in  the  horizontal  and  vertical  components  of  stress 
versus  error  in  the  gage  pattern  radius.  The  results  are 
based  on  the  assumption  that  each  of  the  three  strain 
gages  is  placed  the  same  distance  from  the  center  of  the 
hole.  Fig.  29  essentially  shows  that  the  results  in  terms 
of  stress  are  relatively  sensitive  to  errors  in  the  strain 
gage  pattern  radius  and  that  the  horizontal  components  of 
stress  are  more  sensitive  than  the  vertical  component. 


Consequently  care  must  be  taken  in  the  device  to  measure 
accurately  the  radius  of  the  strain  gage  pattern. 

The  relatively  high  sensitivity  of  the  stress  compo¬ 
nents  to  changes  in  the  gage  pattern  radius  leads  one  to 
question  the  sensitivity  of  the  stresses  to  an  error  in 
the  placement  of  one  of  the  gages  in  the  pattern.  The 
vertical  gage  was  chosen  for  study  since  it  probably 
represents  the  worst  condition  if  not  accurately  placed. 
Again  the  same  nominal  stress  field  was  chosen  and  the 
position  of  the  vertical  gage  was  changed  while  the  two 
gages  at  +  120  degrees  were  maintained  at  a  constant 
radius.  Fig. 30  shows  the  per  cent  error  in  the  vertical 
and  horizontal  components  of  stress  that  would  result  for 
various  errors  in  the  placement  of  the  vertical  gage.  It 
is  noted  that  the  vertical  component  of  stress  is  much  more 
sensitive  to  errors  in  th  -  placement  of  the  vertical  gage 
than  the  horizontal  compoi ents  of  stress  are.  Again  we 
note  that  great  attention  should  be  paid  to  the  positioning 
of  the  strain  gages  in  the  pattern  and  the  exact  measure¬ 
ment  of  the  pattern  prior  to  testing. 

Another  potential  source  of  error  in  the  determination 
of  stress  components  from  relieved  strains  involves  the 
accuracy  of  the  material  properties  determined  from  core. 
Since  the  equations  show  that  the  components  of  stress 


may  be  made  dimensionless  by  dividing  by  Young's  modulus, 
any  error  in  the  modulus  will  be  reflected  direct  y  in 
the  stress  components.  The  situation  is  not  so  straight 
forward  when  considering  the  effect  of  errors  in  Poisson's 
ratio.  The  previously  derived  equations  may  be  used  to 
study  the  effect  of  changes  in  Poisson's  ratio  on  che 
components  of  stress  starting  again  with  a  vertical  stress 
of  1000  psi  and  a  horizontal  stress  of  250  psi  which 
implies  in  the  perfectly  elastic  situation  that  Poisson's 
ratio  is  0.2.  Fig  .31  shows  the  error  in  stress  components 
versus  error  in  Poisson's  ratio.  Vertical  stress  is  less 
sensitive  to  the  error  than  horizontal  stress  although 
neither  is  very  sensitive.  For  example,  an  error  of  25% 
in  Poisson's  ratio  leads  to  an  error  of  2.5%  in  the 
vertical  stress  and  approximately  7%  in  the  horizontal 
stress . 

CONCLUSIONS 

An  approximate  analysis  of  strain  relief  data  from  a 
deephole  trepanning  device  is  possible.  This  approximate 
analysis  is  based  on  a  modified  version  of  the  Soete  method 
in  conjunction  with  the  Kirsch  solution.  The  solution 
presented  is  restricted  to  the  case  where  the  vertical 
direction  is  a  principal  direction.  An  extension  of  the 
approach  presented  here  to  a  more  general  3D  stress  field 
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may  be  possible  by  replacing  the  Kirsch  solution  with  the 
solution  presented  by  Hiramatsu  and  Oka  (1962).  Limited 
experimental  data  tends  to  confirm  the  validity  of  the 
basic  approach  and  agrees  reasonably  well  with  the  analysis 
in  view  of  the  potential  errors  in  material  properties  used 
as  well  as  experimental  errors  in  the  measured  strains. 

Solutions  developed  may  be  used  to  provide  input  into 
the  design  of  the  device.  It  has  been  shown  that  the  ratio 
of  strain  gage  grid  length  to  trepanning  hole  diameter 
3hould  be  kept  as  small  as  possible,  2)  positioning  of 
the  strain  gages  accurately  in  the  pattern  is  essential 
and  3)  computed  stresses  are  directly  sensitive  to  errors 
in  Young's  modulus  of  elasticity  but  much  less  sensitive 
to  errors  in  Poisson's  ratio. 
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PROGRAM  SAP 

I  IDENTIFICATION 

SAP  -  Structural  Analysis  Program 

Programmed  by  E.  Wilson,  L.  Jones,  H.  Dovey  and 
T.  Hsueh  (1970) 

Program  modified  for  use  on  CDC-3400  at  the  South 
Dakota  School  of  Mines  and  Technology  by  M.  Sateesha. 

II  PURPOSE 

The  purpose  of  this  computer  program  is  to  determine 
nodal  displacements  and  element  stress  resultants  of  any 
three  dimensional  solid  subjected  to  general  loading  by 
performing  linear  elastic  analysis.  The  definitions  for  the 
input  and  output  data  are  given  below: 

III  INPUT  DATA 

For  each  three-dimensional  structure  to  be  analyzed,  a 
group  of  punch  cards  is  required  in  this  sequence. 

A.  TITLE  CARD  (12A6) 

Columns  1-72  Alphanumeric  title  for  problem  identifi¬ 
cation 

B.  CONTROL  CARD  (315) 

Columns  1-5  Number  of  nodal  points 
6-10  Number  of  element  types 
11-15  Number  of  load  cases 
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ns 

C.  NODAL  POINT  DATA  CARDS  ( 715,  3F10 . 0,  15 ) 

One  card  per  nodal  point.  Nodal  coordinate  cards  need 
not  be  in  nodal  order  sequence.  If  cards  are  omitted  the 
nodal  data  for  a  series  of  nodes  is  generated. 

Columns  1-5  Identification  -  Nodal  number 

6-10  B.C.  Code  for  displacement  in  X-direction 

11-15  B.C.  Code  for  displacement  in  Y-direction 

16-20  B.C.  Code  for  displacement  in  Z-direction 

21-25  B.C.  Code  for  rotation  about  X-axis 

26-30  B.C.  Code  for  rotation  about  Y-axis 

31-35  B.C.  Code  for  rotation  about  Z-axis 

36-45  X-ordmate 
46-55  Y-ordinate 
56-65  Z-ordinate 
66-70  Mesh  generation  parameter 
A  boundary  condition  code  of  zero  or  blank  indicates 
that  the  joint  is  free  to  move  in  that  direction  and  ^oads 
may  be  applied.  A  boundary  condition  code  of  one  indicates 
that  the  joint  is  fixed  in  that  direction. 

If  a  particular  degree  of  freedom  is  fixed  for  a  series 
of  cards  this  may  be  indicated  by  a  boundary  condition  code 
of  -1  on  the  first  card  in  the  series  and  +1  on  the  last 


card  in  the  series. 
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D.  ELEMENT  DATA  CONTROL  CARD  (415) 

Columns  1-5  The  number  5 

6-10  Total  number  of  elements 
11-15  Number  of  different  materials 
16-20  Number  of  element  distributed  load  sets 

E.  MATERIAL  PROPERTY  CARDS  (I5,4F10.0) 

Columns  1-5  Material  identification  number 
6-15  Modulus  of  elasticity,  E 
16-25  Poisson's  Ratio, 

26-35  Weight  density  of  material 
36-45  Coefficient  of  thermal  expansion 

F.  ELEMENT  DISTRIBUTED  LOAD  SET  CARDS  ( 21 5 , 2F10 . 0, 15 ) 

Columns  1-5  Load  set  identification  number 

6-10  Load  type:  1  for  constant  surface 

pressure,  2  for  hydrostatic  pressure 
11-20  Pressure  for  load  type  1  or  specific 
gravity  for  load  type  2 
21-30  Reference  water  level 
31-35  Element  face  which  pressure  acts  upon 

G.  REFERENCE  TEMPERATURE  (2F10.0)  CARD 

Columns  1-10  Stress  free  temperature 

11-20  Acceleration  due  to  gravity 

H.  ELEMENT  LOAD  CASE  FACTOR  CARDS  (5  cards  of  4F10.0) 

Pressure  and  thermal  load  factors  on  the  element  load 


> 
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cases  are  scaling  factors  in  order  to  provide  flexibility 
in  modifying  applied  loads. 


Card 

1: 

Columns  1-10 

Pressure  load 

A 

11-20 

factors  for 

B 

21-30 

element  load 

C 

31-40 

cases 

D 

Card 

2: 

Columns  1-10 

Thermal  load 

A 

11-20 

factors  for 

B 

21-30 

element  load 

C 

31-40 

cases 

D 

Card 

3: 

Columns  1-10 

Percentage  of  gravity 

A 

11-20 

acting  in  +X 

B 

21-30 

direction  in  element 

C 

31-40 

load  case 

D 

Card 

4: 

Columns  1-10 

Percentage  of  gravity 

A 

11-20 

acting  in  +Y 

B 

21-30 

direction  in  clement 

C 

31-40 

load  case 

D 

Card 

5 : 

Columns  1-10 

Percentage  of  gravity 

A 

11-20 

acting  in  +Z 

B 

21-30 

direction  in  element 

C 

G.  ELEMENT  CARDS  ( 1215 , 412, 211, F10 . 2) 
Columns  1-5  Element  number 


C-5 


6-10 

Global 

1 

11-15 

node 

2 

16-20 

point 

3 

21-25 

numbers 

4 

26-30 

corresponding 

5 

31-35 

to 

6 

36-40 

element 

7 

41-45 

nodes  - 

8 

46-50 

Integration  order 

51-55 

Material  number 

56-60 

Generation  Parameter 

61-62 

Distributed  load  set 

A 

63-64 

number  for 

B 

65-66 

element  load  cases 

C 

67-68 

(zero  implies  no  load) 

D 

69-70 

Face  numbers  for  stress  output 

71-80 

Element  temperature 

Note:  1.  Element  cards  must  be  in  ascending  order 

2.  Computation  time  increases  with  the  cube  of  the 
integration  order.  The  order  is  2  for  rectangular 
elements  and  3  for  skewed  elements 
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3.  Element  faces  are  numbered  as  follov/s 
Face  1  Corresponds  to  + £  direction 

2  Corresponds  to  -  £  direction 

3  Corresponds  to  +  n  direction 

4  Corresponds  to  - n  direction 

5  Corresponds  to  +  q  direction 

6  Corresponds  to  -  5  direction 

0  Corresponds  to  the  center  of  the  element 

4.  Sign  convention:  A  positive  distributed  load  acts 
in  the  positive  (local)  axis  direction  associated 
with  each  face. 

H.  CONCENTRATED  LOAD  DATA  CARDS  (2I5.6F10.0) 

One  card  per  load  case  for  each  node  which  has  non¬ 
zero  concentrated  loads  or  moments  applied.  The  cards  must 
be  in  nodal  number  sequence 

Columns  1-5  Nodal  number 

6-10  Load  condition  number 
11-20  Load  in  X-direction 
21-30  Load  in  Y-direction 
31-40  Load  in  Z-direction 
41-50  Moment  about  X-axis 
51-60  Moment  about  Y-axis 
61-70  Moment  about  Z-axis 


J.  BLANK  CARD 


) 
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) 


The  above  sequence  of  cards  must  be  terminated  with 
one  blank  card. 

K.  ELEMENT  LOAD  MULTIPLIER  CARDS  (4F10.0) 

One  card  must  be  supplied  for  each  load  condition 
which  contains  the  following  information 

Columns  1-10  Multiplier  for  element  load  A 
11-20  Multiplier  for  element  load  B 
21-30  Multiplier  for  element  load  C 
31-40  Multiplier  for  element  load  D 
These  cards  must  be  in  load  order  sequence. 

IV  OUTPUT  INFORMATION 

The  program  prints  the  following  output 

A.  Input 

B.  Nodal  displacements 

C.  Element  Stresses. 


PROGRAM  SOLID 


I  IDENTIFICATION 

SOLID  -  A  modified  version  of  SAP 

Program  developed  for  use  on  CDC-3400  at  the  South 
Dakota  School  of  Mines  and  Technology  by  M.  Sateesha. 

II  PURPOSE 

The  purpose  of  this  program  is  to  generate  the  element 
stiffness  and  stress  matrices  for  eight  node  hexahedral 
three-dimensional  solid  elements.  The  definitions  for  the 
input  and  output  data  are  given  below. 

III  INPUT  DATA 

Same  as  in  sections  A-G  as  described  under  Program  SAP. 

IV  OUTPUT  INFORMATION 

A.  Element  stiffness  and  pointer  matrices  on  tape 
unit  2 . 

B.  Element  stress  matrices  on  tape  unit  8. 
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PROGRAM  ASEMBLE 

I  IDENTIFICATION 

ASEMBLE  -  Assemblage  of  element  stiffness  matrices 

Programmed  by  M .  K.  Sateesha 

Department  of  Geological  Engineering 
South  Dakota  School  of  Mines  and  Tech¬ 
nology,  Rapid  City 

II  PURPOSE 

The  purpose  of  this  computer  program  is  to  assemble 

the  element  stiffness  matrices  into  a  global  stiffness 

matrix. 

III  INPUT  DATA 

A.  DISK  CONTROL  CARD  (215) 

Columns  1-5  Maximum  number  of  records  on  random 
access  disk 

6-10  Maximum  number  of  rows  that  can  be 
handled  at  a  time. 

B.  INPUT  TAPE  ON  UNIT  2 

This  tape  is  obtained  as  output  from  program  SOLID. 

It  contains  the  element  stiffness  and  pointer  matrices. 

C.  CONTROL  CARD  (215) 

Columns  1-5  Number  of  equations 

6-10  Number  of  nodes  at  which  load  boundary 


conditions  are  specified. 


D.  LOAD  CONDITION  CARDS  (4I5.3F10.4) 


Columns  1-5  Nodal  number 

6-10  B.C.  Code  for  displacement  in  +X  direction 

11-15  B.C.  Code  for  displacement  in  +Y  direction 

16-20  B.C.  Code  for  displacement  in  +  Z  direction 

21-30  Load  in  X-direction 
31-40  Load  in  Y-direction 
41-50  Load  in  Z-direction 


IV  OUTPUT  INFORMATION 

The  program  prints  the  input  and  writes  the  assembled 
stiffness  matrix  (global),  load  vector  and  pointer  matrix 
on  tape  unit  4. 
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PROGRAM  GSOR 

I  IDENTIFICATION 

GSOR  -  Gauss-Seidel  over  Relaxation  Technique 
Programmed  by  M.  K.  Sateesha 

Department  of  Geological  Engineering 
South  Dakota  School  of  Mines  and 
Technology,  Rapid  City 

II  PURPOSE 

The  purpose  of  this  computer  program  is  to  solve  a 
set  of  linear  simultaneous  equations  utilizing  the  well 
known  Gauss-Seidel  iteration.  The  solution  is  refined  at 
the  end  of  every  cycle  by  using  an  over-relaxation  factor. 

The  program  is  specifically  suited  for  banded  equations 
since  the  zero  coefficients  are  not  stored  and  the  corres¬ 
ponding  operations  skipped  thus  saving  both  storage  and  time. 
Hi  INPUT  DATA 

A.  INPUT  PARAMETER  CARD  <2I5,F10.0) 

Columns  1-5  Number  of  equations  to  be  solved 
6-10  Maximum  number  of  iterations 
11-20  Tolerance  allowed 

B.  INITIAL  DISPLACEMENT  CARDS  (4E20.8)  -  OPTIONAL 

If  these  cards  are  not  provided,  the  initial  values 


are  all  assumed  to  be  zeros. 
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IV  OUTPUT  INFORMATION 

The  program  prints  the  following  output: 

A.  The  starting  displacement  values. 

B.  Whether  the  solution  converged  or  not. 

C.  The  number  of  iterations  required  for  converging  to 
the  final  solution. 

D.  The  final  solution. 


PROGRAM  DISTRES 


I  IDENTIFICATION 

DISTRES  -  Determination  of  element  stresses  and  strains 
Programmed  by  Malalur  K.  Sateesha 

Department  of  Geolog ical  Engineering 
South  Dakota  School  of  Mines  and  Tech¬ 
nology,  Rapid  City 

II  PURPOSE 

The  purpose  of  this  program  is  to  compute  the  element 
stresses  and  strains  from  the  nodal  displacements. 

III  INPUT  DATA 

A.  TITLE  CARD  (12A6) 

Columns  1-72  Alphanumeric  title  for  problem  identifi¬ 
cation 

B.  CONTROL  CARD  (215) 

Columns  1-5  Number  of  nodal  points 
6-10  Number  of  element  types 

C.  MATERIAL  PROPERTY  CARE  (2F10.0) 

Columns  1-10  Young's  Modulus 
11-20  Poisson's  Ratio 

D.  NODAL  DISPLACEMENT  CARDS  (4E20.8) 

These  cards  are  obtained  as  output  from  Program  GSOR. 

E.  TAPE  UNIT  1  -  OUTPUT  TAPE  UNIT  8  FROM  PROGRAM  SOLID 


Ill  OUTPUT  INFORMATION 

1 .  Input  Data 

2.  Nodal  Displacements 

3.  Element  Stress  Resultants 

4.  Element  Strain  Resultants 
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PROGRAM  PRIN 

I  IDENTIFICATION 

PRIN  -  Principal  Stress  Determination  Program 

II  PURPOSE 

The  purpose  of  this  computer  program  is  to  determine 
the  principal  stresses  and  their  directional  cosines  for  any 
given  stress  tensor. 


Ill  INPUT  DATA  (3  cards  of  3F20.8) 

The  stress  tensor  is  read  rowwise  and  stored  as 


x 


xy 


xz 


XX 


xz 


zx  zy 

IV  OUTPUT  INFORMATION 


The  principal  stresses  are  obtained  as 


0  l 

0 

0 

0 

02 

0 

0 

0 

03 

and  the 

directional  cosines  as 

Vl-X 

V2-X 

V3-X 

Vl-Y 

V2-Y 

V3-Y 

where 

Vl-Z 

V2-Z 

V3-Z 

VI, 

V2,  V3 

=  Unit 

vectors 

o3  respectively, 
and 


X,  Y,  Z  =  components  of  the  vectors  along  the  Global 


coordinate  directions 


LIST  OF  REFERENCES 


I3<> 


ARGYRIS,  J.H.,  1955,  Energy  Theorems  and  Structural 

Analysis,  Aircraft  Engineering,  v.  27,  p.  125-154. 

CI/5UGH,  R.W.,  1965,  The  Finite  Element  Method  in  Struc¬ 

tural  Mechanics,  Stress  Analysis  (Ed.  Zienziewicz 
and  Hollister) .  John  Wiley  and  Sons,  New  York, 

Chap.  7,  p.  85-119. 

DANIELLS,  P.#  1971,  Experimental  Measurement  of  Residual 

Strain  in  Rocks,  M.S.  Thesis,  S.D.  School  of  Mines 
and  Technology,  Rapid  City. 

EMERY,  C.L.,  1964,  "Strain  Energy  in  Rocks,"  State  Of  Stress 

in  Earth's  Crust,  American  Elsevier,  New  York, 
pp.  234-279. 

FELIPPA,  C.A.  and  CLOUGH,  R.W.,  1970,  The  Finite  Element 

Method  in  Solid  Mechanics,  Numerical  Solution  of 
Field  Problems  in  Continuum  Physics,  SIAM-AMS  Proc., 
v.  II,  Amor.  Math.  Soc . ,  Providence,  Rhode  Island. 

FRIEDMAN,  M.,  1968,  "X-ray  Analysis  of  Residual  Elastic 

Strains  in  Quartzose  Rocks":  Tenth  Sym .  on  Rock  Mech . , 
(in  press),  Austin,  Texas. 

GOLDICH,  S.S.  et.al.,  1961,  The  Precambrian  Geology  and 
Geochronology  of  Minnesota,  Minnesota  Geological 
Survey  Bull.,  41,  Minneapolis,  193  p. 

HAFNER,  W. ,  1951,  Stress  Distribution  and  Faulting, 

Geol .  Soc.  America  Bull.,  v.  62,  p.  373-398. 

HAIMSON,  B.,  1973,  Real  Stresses  Around  Boreholes, 

Proceedings  of  the  Sixth  Conference  on  Drilling  and 
Rock  Mechanics,  The  University  of  Texas,  Austin,  Texas. 

HAIMSON,  B.,  LaCOMB,  J.,  JONES,  A.H.,  and  GREEN,  S.J.,  1974, 
Crustal  Stresses  At  The  Nevada  Test  Site,  Transaction 
of  The  A.G.U.,  vol.  55,  no.  4,  p.  425. 


HALL,  1899,  The  Gneisses,  Gabbros,  Schist  and  Associated 
Rocks  of  Southwestern  Minnesota,  USGS  Bull,  157. 


J 


131 


HEALY,  J.H.,  ROLLER,  J.C.,  GORTON,  Q.,  and  LAMSON,  R. 

1974,  In  Situ  Stress  Measurements  in  Fractured  Rocks, 
Transactions  of  the  A.G.U.,  vol .  55,  no.  4,  p.  425 

HIRAMATSU,  Y.  and  OKA,  Y.,  1962,  Stress  Around  a  Shaft 

or  Level  Excavated  in  Ground  with  a  Three  Dimensional 
Stress  State,  Faculty  of  Engineering  Memoirs,  Kyoto 
University,  Japan,  24. 

HOSKINS,  E.,  1967,  An  Investigation  of  Strain  Rosette  Relief 
Methods  of  Measuring  Rock  Stress,  Int .  J .  Rock  Mech . 

Min.  Sci.,  v.  4,  p.  155-164. 

HOSKINS,  E.,  1968,  A  Stress  Measurement  Technique  fcr  Deep 

Holes,  Trans.  AGU.v.  49,  no.  4,  p.  756. 

HOWARD,  J.H.,  1966,  Structural  Development  of  the  Williams 

Range  Thrust,  Colorado,  Geol .  Soc ,  America  Bull.,  v.  77 
p.  1247-1264. 

HOWLAND,  R.C.J.,  1930  On  the  Stresses  in  The  Neighborhood 

of  a  Circular  Hole  on  a  Strip  Under  Tension,  Phil.  Trans., 
Roy.  Soc.  (London),  Series  A.  v.  229,  p.  49-86. 

KEHLE,  R.O.,  1964,  Determination  of  Tectonic  Stresses 

Through  Analysis  of  Hydraulic  Well  Fracturing,  J.  Geoph . 
Res.,  v.  69.  p.  259-273. 

KIRSCH.,  G.  1898,  "Die  Theorie  der  Elastizitat  und  die 
Bedurfnisse  der  Festigkeitslehre „  "  Zeitschr ift 
Verein  Deutscher  Ingenieure:  vol.  42,  no.  29  pp. 

797-807. 

McCLINTOCK,  F.A.  and  ARGON,  A.S.,  1966  Mechanical  Behavior 
of  Materials,  Addison-Wesley  Pub.  Co.,  Reading,  Pa. 

PRZEMIENIECKI ,  J.S.,  1968,  Theory  of  Matrix  Structural 

Analysis,  McGraw  Hill  Book  Co.,  New  York,  468  p. 

SANFORD,  A . R . ,  1959,  Analytical  and  Experimental  Study 

of  Simple  Geologic  Structures,  Geol.  Soc.  America 
Bull . , v .  70,  p.  19-52. 

SATEESHA,  M.K.,  1974,  "Analysis  of  Data  from  a  Deep  Hole 

Stress  Measurement  Device,"  Ph.D.  thesis,  S.D.  School 
of  Mines  and  Technology. 


132 


SOETE,  W.  and  VANCROMB RUGGE,  R.,  1950,  "An  Industrial  Method 
for  the  Determination  of  Residual  Stresses,"  Proc .  - 
Society  for  Experimental  Stress  Analysis,  vol.  8,  no.  1, 
pp.  17-28. 

THEIL,  G.A.  and  DUTTON,  C.E.,  1935,  The  Architectural, 

Structural  and  Monumental  Stones  of  Minnesota, 

Minnesota  Geol.  Survey,  Bulletin  No.  25,  Minneapolis, 

p.  160. 

VARNES,  D.J.,  1969,  "Model  for  Simulation  of  Residual  Stress 
in  Rock,"  Rock  Mechanics-Theory  and  Practice.  Proc. 
Eleventh  Sym.  on  Rock  Mech.,  AIME,  New  York,  pp.  415-426. 

WHITE,  JON,  1973,  Effects  of  Small  Temperature  Changes 
on  Measured  Apparent  Residual  Strain  in  Selected 
Rock  Cylinders,  M.S.  Thesis,  South  Dakota  School  of 
Mines  and  Technology,  Rapid  City,  South  Dakota. 

ZIENKIEWICZ,  O.,  1971,  The  Finite  Element  Methods  in  Engineering 
Science.  McGraw  Hill  Book  Co.,  London,  p.  521. 


